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ABSTRACT 


This final report summarizes the work performed under Exhibit A 
of NASA Contract NAS9-12183. It includes the simulation developments 
for use in dynamics and control analysis during boost from liftoff 
to orbit insertion. It also includes wind response studies of the 
NR-GD 161B/B9T delta wing booster/delta wing orbi ter configuration, 
the MSC 036B/280 inch solid rocket motor configuration, the MSC 
040A/L0X-propane liquid injection TVC configuration, the MSC. 040C/ 
dual solid rocket motor configuration-, and .the MSC 049/solid rocket 
motor configuration. All of the latest math models (rigid and flexible 
body) developed for the MSC/GD Space Shuttle Functional Simulator, 
are included in this report. 
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1.0 


SUMMARY 


This final report is submitted to the NASA MSC Guidance and Control Division 
as required by MSC DRD SE-272T of Exhibit A of Contract NAS9-12183. This 
report covers the work performed from contract go-ahead to contract completion 
on 15 September 1972. Under the terms of the contract, analysis programs were 
to be updated (Task I), rigid body dynamic and control analyses were to be 
performed (Task II), and flexible body stability and control analyses were to 
be performed (Task III). A brief summary of these studies, which have been 
reported to MSC/GCD in 49 interim reports (listed in the references) is given 
in the following paragraphs. , 

Task I Summary - Update Of Analysis Programs 

Initially, mathematical models were developed for the MSC Space Shuttle 
Functional Simulator to provide a launch phase boost to staging simulation 
equivalent to the Boeing-developed Rigid Body Boost Simulation (RIBBS). 

These models were programmed by MSC C&AD into the Space Shuttle Functional 
Simulator (SSFS) developed by the Guidance and Control Division. The SSFS 
first stage simulation was checked out jointly by MSC and Boeing personnel. 

In order to provide analysis capability from liftoff through orbit insertion, 
several additional models not already in RIBBS were developed for SSFS. These 
models were: targeting to orbit; active guidance; a three degree of freedom 

upper stage to orbit; and calculation of orbital parameters. During the 
course of the model development, programming, and checkout of SSFS, numerous 
modifications were made to the RIBBS program as a result of updated analysis 
requirements in Task II. SSFS was then updated to meet these additional 
analysis capability requirements by adding: trajectory shaping and control 

gain calculations; baseline control system, multistage guidance; and liftoff 
state vector calculations. 

Point time stability analysis equations were also developed during Task I 
and programmed into the M DELTA Program. 

All of the math models which were developed for the SSFS have been updated 
and are included in this final report. 
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Task II Summary - Rigid Body Dynamics and Control Analysis - 

The NR/GD 161B/B9T delta wing booster/delta wing orbiter configuration was 

extensively analyzed during the first half of the contract. Brief studies 

were performed on the 036B/280 inch solid rocket motor series burn configuration 

and the 040A/L0X propane series burn configuration using liquid injection thrust 

vector control. In February 1972 studies. were initiated on parallel burn 

dual solid rocket motor (SRM) configurations. 

Liftoff studies of the 040C/SRM configuration revealed that 0.25° SRM mis- 
alignments and ground winds were satisfactorily controllable. However, 0.5° 
misalignments caused thrust torques exceeding the TVC roll control authority. 
Orbiter engine out analysis demonstrated that the SRMS must be oriented to 
thrust nearly through the liftoff eg or the vehicle would be uncontrollable - 
in event of an orbiter engine failure. Inflight analysis of the 040C/SRM 
showed that the available roll torque was insufficient to maintain a fixed 
roll attitude in crosswinds, so studies of active roll control to minimize 
sideslip were conducted. 

The configuration was changed to the MSC049 Model and control system develop- 
ment was continued. The resulting control system is based on a gravity turn 
trajectory so the vehicle can roll and keep trajectory dispersions to a 
minimum. The roll control problem is minimized by sensing sideslip and 
actively rolling its tail into the wind. Results of several hundred simulated 
flights to orbit indicated that this control system can handle any wind 
conditions as Tong as SRM misalignments are constrained within 0.25 degrees. 
Aerodynamic loads were automatically minimized without load relief techniques, 
and weight to orbit was maintained within a few thousand pounds of the nominal 
no-wind value for all cases. 

Because of the 0.25 degree alignment constraint studies of trim gimballing 
of the SRMS were also performed. It was determined that gimballing the SRMS 
+ two degrees at a rate limited to 0.3 degrees/second not only permitted 
misalignments greater than one degree but also improved wind response of 
the vehicle. 
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Task III Summary - Flexible Body Stability and Control Analysis 
Flexible body models were developed and programmed into the space shuttle 
functional simulator so that flexible body analysis can be initiated as 
soon as bending data becomes available. The flexible body program uses 
a generalized modal approach to bending which represents the elastic response 
by standard normal modal equations with viscous damping. Included are models 
for distributed aerodynamic forces and moments, and thrust forces and moments 
to account for bending effects as well as the tail wags dog contribution to 
bending. The vibration model sums all the forces and moments acting on each 
of the equivalent mass points and for a given mode numerically integrates 
the sum with a second order differential equation in modal displacement. 
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2.0 Task I - Update of Analysis Programs 

The three major programs which were included as part of the updating procedure 
of this task were the Space Shuttle Functional Simulator (SSFS), Rigid Body 
Boost Simulation (RIBBS), and Matrix Differential Equation Linear Transform 
Analysis (MDELTA). SSFS was updated with mathematical models to provide a 
launch phase simulation for the MSC Guidance and Control Division. All of 
the models which were developed during this contract period have been updated 
to the current configuration and are. presented in Section 2.1. The Boeing 
RIBBS program was used for analysis of shuttle during the development of 
SSFS. RIBBS was also updated with a new guidance model and other modifications 
as explained in Section 2.2. The Boeing MDELTA program was modified for use 
in point time stability analysis. This is mentioned in Section 2.3. Also 
included in Section 2.4 is a discussion of the development of the boost 
Linear Tangent Guidance (LTG) equations. 
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2.1 


Space Shuttle Functional Simulator (SSFS) Math Models 

Rigid body mathematical models were prepared and delivered to MSC Guidance 
and Control Division to provide the Space Shuttle Functional Simulator (SSFS) 
with the same capability as the Boeing Rigid Body Boost Simulation (RIBBS). 
These models which were originally presented in References 1 and 2 include: 
boost polynomial guidance; mass properties thrust; thrust vector control ; 
rigid body 6 DOF equations of motion; aerodynamics; RCS; and atmosphere 
models. These models were programmed by Lockheed Electronics Corporation 
(LEC) and jointly checked out by LEC and Boeing by making comparison runs 
on RIBBS and SSFS. References 3 and 4 describe the comparison results. 

Slight differences were observed due to some fundamental program structure 
differences and updated gravity and atmosphere models in SSFS. The final 
checkout used design shear and gust winds to verify the short period dynamics 
and included detailed examination of the computer programming to assure that 
the differences that remained after debugging were not the result of coding 
errors. Additional confidence was gained by duplicating some of the more 
complex transformations and equation sets on the Boeing computers (using the 
Boeing computer service remote terminals) and verifying them using known 
input-output data sets. 

Severaladditional new models were also developed and delivered to MSC. 

These include orbit insertion targeting, active guidance, three degree of 
freedom upper stage to orbit equations of motion, orbital parameters, 
initial state vector, hinge moment calculations, control gains calculations, 
and trajectory shaping. These models were detailed in References 5-11. 

These models were also programmed by LEC and comparisons of results with 
existing simulation data was favorable. The control gains and trajectory 
shaping equations were verified, as shown in References 12 and 13, by 
comparison with similar equations programmed into the RIBBS program. 

The following subsections give a complete description of the current math 
models developed during this contract and incorporated in the Space Shuttle 
Function Simulator. The baseline control system of Reference 14 that was 
accepted for inclusion in the MSC Guidance and Control Equations Document 
"MSC-04217 Revision B" has since been updated to include sideslip feedback and 
gimballed solid rocket motors as presented in paragraph 2.1.5. 
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2.1.1 AERO (Aerodynamics) 

2. 1.1.1 DESCRIPTION 

This model calculates the latitude and longitude of the vehicle 
position, the contribution to velocity due to wind speed and 
direction, flight path angle, Mach number, dynamic pressure, 
angle of attack and angle of sideslip. The vehicle forces and 
moments due to the air and control surface deflections are cal- 
culated and summed. 
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2.1.1 .AERO (Continued) 
2. 1.1. 2 MATH MODEL 
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2 . 1.1 


AERO (Continued) 


V w = table lookup ~f (altitude) 
AZ W = table lookup ~f (altitude) 




Obtain from ATM0S 

1) p = f (altitude) 

2) a = f (altitude) 

3) P = f (altitude) 
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2.1.1 AERO (Continued) 


M = Vg/a 



The following aerodynamic coefficients are looked up in tables 

C 7 = f (M) 
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2 . 1.1 

Where: 


AERO (Continued) 




V EARTH X 

V EARTH y 






LV 



= vehicle position in inertial polar-equatorial coordinates 

= vehicle position in inertial plumbline coordinates 

= transformation matrix from inertial polar-equatorial to 
plumbline coordinates 

= latitude of present position of vehicle 

= East longitude of present position of vehicle corrected 
for earth's rotation 

= time of launch (from epoch) i 

= radians to degrees conversion constant 

= angular rate of earth 

= elapsed time from liftoff 

= distance from the center of the earth to the yehicle 

= components of earth's velocity in inertial polar- 
equatorial coordinates 

= components of vehicle relative velocity in plumbline coordinates 


= components of vehicle velocity in plumbline coordinates 
= total vehicle relative velocity in plumbline coordinates 

= components of relative velocity in local vertical coordinates 

= transformation matrix from local vertical to inertial polar- 
equatorial coordinates 


Y = vehicle flight path angle with respect to local horizontal 
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2 . 1.1 


u H 



M 


a 

3 


P 

a 

P 

M 



b 


AERO (Continued) 

= horizontal wind speed in local vertical coordinates 
= wind azimuth (North = 0°) 

P = components of wind velocity in plumbline coordinates 

3 

B = vehicle velocity with respect to air in body coordinates 

= transformation matrix from body to plumbline coordinates 

= total vehicle velocity with respect to air in body coordinates 

= vehicle angle of attack 
= vehicle sideslip angle 
= local air mass density 
= local speed of sound 
= local air pressure 
= Mach number 
= dynamic pressure 

= components of aerodynamic force in body coordinates 

= vehicle aerodynamic reference area 
= vehicle mean aerodynamic chord 
= vehicle reference span 
= aileron deflection 
= elevator deflection 
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2.1.1 


AERO (Continued) 


6 r = rudder deflection 

P = vehicle roll rate 

Q = vehicle pitch rate 

R = vehicle yaw rate 

5 

X Y Aerodynamic moments about the X, Y, and Z body axes, 
respectively 
M Z 
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IU 


2 . 1.1 


AERO (Continued) 


2. 1.1. 3 INPUT/OUTPUT ++ 

Input from routines: 



vehicle position in inertial plumbline 
coordinates from EOM 



vehicle velocity in inertial plumbline 
coordinates from EOM 



wind velocity from tables 



wind azimuth from tables 


p » a » p 


current air density, speed of sound and air 
pressure from ATMOS 




current location of vehicle center of gravity 
from MASPRO 


t 


elapsed time from liftoff from flight sequencer 


distance from center of the earth to the vehicle 
from E0M 


W 6 R 


aerodynamic control surface deflections from 
flight software commands 


P, Q, R 


vehicle roll, pitch and yaw rates from E0M 


All aerodynamic coefficients are input from tables. 


Input from cards for initialization: 



b 

X A • y A ’ Z A 
Mr Mr Mr 


time of launch (from epoch) 

radians to degrees conversion constant 

angular rate of earth 

vehicle aerodynamic reference area 

vehicle mean aerodynamic chord 

vehicle reference span 

aerodynamic reference location in body 
coordi nates 


1.4 



2.1.1 


AERO (Continued) 


Output to routines: 


tt X h Y m Z 

^A » ) ^A 

M X rt Y R 1 


Output to printer: 
4» 

v x 

X LV 


Y 

V W 

AZ 


W 


'B 


a 

B 

P 

a 

P 

M 

q 


V V 6 a 


current air pressure to THRUST 
components of aerodynamic forces to the E0M 

moments due to aerodynamic forces to E0M 

latitude of vehicle's position 
longitude of vehicle's position 
rate of cl imb 

flight path angle 
wind speed 
wind azimuth 

vehicle velocity with respect to air 

angle of attack 

angle of sideslip 

local air mass density 

local speed of sound 

local air pressure 

Mach number 

dynamic pressure 

aerodynamic control surface deflections 
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2.1.2 ATMOS (Atmosphere Model) 

2. 1.2.1 PROGRAM DESCRIPTION 

This program calculates the speed of sound, pressure and air density 
from an altitude input. 

2. 1.2.2 MATH MODEL 

Use the Cape Kennedy Reference Atmosphere (TM-X-53872, PARAGRAPH 14.7 
MSFC "COMPUTER SUBROUTINE PRA-63") as specified for SSV design studies 

2. 1.2.3 INPUT/OUTPUT 

The altitude above the mean earth surface must be supplied to the 
model which returns the speed of sound, pressure, and atmospheric 
density. 
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2.1.3 ATTUDE (Boost Polynomial Guidance) 

2. 1.3.1 PROGRAM DESCRIPTION 


This program is used in the flight software to guide the vehicle 
during atmospheric flight or until active guidance is commanded to 
take over the guidance function. The yaw and roll commands are 
set to zero and the pitch command is taken from a pitch table 
look-up. The pitch table is derived by shaping the desired tra- 
jectory (using T SHApE ), and using the pitch angles from this 
trajectory as the commands for the Boost Polynomial Guidance. 


2. 1.3.2 MATH MODEL ; 

* . • ■ 

The pitch commands (TA) and their corresponding time (T) are taken 
from a desired trajectory and entered into a table of commands and 
times. The following logic is then used to derive the commanded 
vehicle attitude during this part of the flight. 
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2 . 1.3 


ATTUDE (Continued) 



where T = time after launch, seconds 
e c = commanded pitch angle 
= 0 . = commanded yaw angle 
<j> c = 0. = commanded roll angle 

The coordinate systems used in this model are the platform systems and the 
body systems which will be described in section IV. 
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2.1.3 ATTUDE (Continued) 

2. 1.3. 3 INPUT/OUTPUT 

This model requires time after launch as an input and it~cal culates the 
roll, pitch, and yaw attitude angles. The program requires no software 

command, calling arguments, or other subroutine. It should be called 
about every two seconds. - - - • " 

2.1. 3. 4 COORDINATE SYSTEMS- •' - ' ' : 

The angles calculated by this model relate the body coordinate system to 
the platform coordinate system, the body system is fixed with respect to 
the vehicle with the X axis forward thru the main propellant tank centerline, 
the Z axis in the engine gimbal pivot plane and pointing down, and the Y axis 
points toward the pilot's right forming the right handed triad. The platform 
system origin is at the earth's center of mass and is fixed in inertial 
space at the time of launch.: The X axis is parallel but opposite in sense to 
the launch pad gravity vector, the Z axis points dcwnrange in the launch plane 
and the Y axis points toward the pilot's right completing the right handed triad. 


; - 

cose 

0 

sine 


[e] r: 

- a 

1 

0 

( 1 ) 


-sine 

0 

cose 



COSi}) 

-sim{) 

0 


M - 

s i nxp 

COSiJ) 

0 

( 2 ) 


0 

0 

1 
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2.1.3 ATTUDE (Continued) 



0 0 
cos 4 > -si n 4 > 

sin<> cos<{> 


T 

Q. 

X 

i 


X B 

y p 

[e] M '■[♦] 

y b 

1 

O- 

M 

1 


Z B 


(3) 



where: X p> Y pJ Z p = platform coordinate vector 

Xg» ’ Z B = body coordinate vector 


2.1.4 


AVEH (6 DOF Equations of Motion) 


2. 1.4.1 DESCRIPTION 

AVEH defines the motions of the center of gravity of the vehicle. 

For convenience it is separated into three parts; 1) translation 
equations, 2) rotation equations, and 3) euler angles. 

These equations should be solved at least once each second during 
powered flight. In the vicinity of environmental discontinuities 
more frequent solution is required; for instance, the vehicle can 
fly completely through a wind gust at maximum dynamic pressure 
within 0.1 second. Other discontinuities include: staging, start 

of closed loop guidance, and engine or actuator failures. As a rule 
of thumb, the integration rate during transients can be % x rotational 

O 

acceleration (in degrees/sec ). 

2. 1.4. 1.1 TRANSLATION EQUATIONS 
DESCRIPTION 

This model defines the linear accelerations of the rigid body. 


k- - - 

EF V 

r 


ZF y 

iFw 

P 

= [B] 

ef y 

B 

1 

M 

-n 

rvi 

“O 

1 


IF 


MATH MODEL 



2.1.4 AVEH (Continued) 


I 

a. 

X 

lZL 



g Y 

p 

= [«] 

\ 

9 Z 

. P L 

- 

1 

t — i 

cn , 

1 


X P = % + ZF x/” *p = / tj * dt + *P 

f p - J^YdtVtp 

z P hi Ut * 2 P 


Y p — g v zFy /m 

Y p T p 


Zp - g 7 + zF 7 /m 
K *P ^P 


o 


o 


0 



ft* ^ + 


f l 2 

>h 

f\\ 


Xdt + Y r 


Zdt + Z r 


o 


o 


o 


IF Y , eF v , zF 7 = sum of forces in the X, Y, Z body axis directions. 
a B y B z B 


eF y , zF v , zF 7 = sum of forces in the X, Y, Z inertial plumbline axis 

Ap Yp L p 

directions. 


[f] = transformation matrix from body to inertial plumbline. 

zF = aero forces + thrust forces + RCS forces + engine deflection 
forces + slosh forces. 



= gravitational acceleration components in inertial polar 
equatorial axis directions. 



= gravitational acceleration components in inertial 
plumbline axis directions. 


[a] = transformation matrix from inertial polar - equatorial to 

inertial plumbline. 


Xp> Yp, Zp = accelerations in inertial plumbline axis directions 
m = total vehicle mass 
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eF y , zF v , iF 7 = sum of forces in the X, Y, Z body axis directions. 
a B t B B 


eF y , zF v , zF 7 = sum of forces in the X, Y, Z inertial plumbline axis 

Ap Tp £p 

directions. 

[B] = transformation matrix from body to inertial plumb line. 

zF = aero forces + thrust forces + RCS forces + engine deflection 
forces + slosh forces. 

g x , g Y , g 7 = gravitational acceleration components in inertial polar- 
A I T I equatorial axis directions. 

g x , gy » 9z = gravitational acceleration components in inertial 
A P T P l P plumbline axis directions: • 

[a] = transformation matrix from inertial polar - equatorial to 

inertial plumbline. 

Xp> Yp, Zp = accelerations in inertial plumbline axis directions 
m = total vehicle mass 
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2.1.4 


AVEH (Continued) 



INPUT/OUTPUT 

The translation equations require as inputs: 

Aerodynamic forces 
Thrust forces 
RCS forces 

Engine deflection forces 

Slosh forces 

[a] and [e] matrices 

Gravitational acceleration components 

Vehicle mass 

Initial conditions on Xp, Yp, Zp, Xp, Y p , Zp 

The outputs from the translation equations are: 

EF Y , ZF V , zF 7 and the inertial plumbline position, velocity 
a B y B ^B 
and acceleration components. 

The translation equations require the presence of subroutines: RCS, 
THRUST, AERO JVC, SLOSH AND GRAVITY, 


COORDINATE SYSTEMS 

Body axes - Orthogonal system with origin at engine gimbal pivot plane - 
X axis positive toward nose of vehicle; Z axis positive "down", and Y 
axis positive toward the right wing. 

Inertial plumbline - Orthogonal system with origin at center of the earth - 
X axis parallel to the launch site gravity vector and positive in the 
direction opposite to gravitational acceleration. 

Transformation matrix from body to inertial plumbline: 
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2.1.4 


AVEH (Continued) 


'11 


'12 


13 


21 

5 22 


'23 


'31 

>32 


'33 


= Cos e Cos \p 

= Sin e Sin 4> -Cos e Sin 41 Cos 4 
= Sin e Cos <j> + Cos 0 Sin 4 Sin 4 
= Sin 4 
= Cos 4 Cos 4 
= -Cos 4/ Sin 4 
= -Sin 0 Cos 4 

= Cos 0 Sin 4 + Sin 0 Sin 4 Cos 4 
= Cos 0 Cos 4 - Sin 0 Sin 4) Sin 4 


(See Euler Angles below 
for calculation of 

Q» *{>><!>) 


Inertial polar - equatorial - orthogonal system with origin at the center 
of the earth - X axis is in equatorial plane, positive through a reference 
meridian at time of liftoff; the reference direction is defined by the 
time of liftoff and the coordinate system used for gravity calculations. 
The Z axis is positive through the North Pole. 

Transformation matrix from polar - equatorial to plumbline 




a ll 

a 12 

a 13 

1 — 1 

p 

1 1 

= 

a 21 

a 22 

a 23 



a 31 

a 32 

a 33 


i 


* 

a, . = Cos A. Cos (4, + W t, ) 

11 L L e L 

a 12 = Cos A l Sin 
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2.1.4 

13 

AVEH (Continued) 
Sin A^ 

" . 

'21 

Sin A1 SinA, Cos (W t, + 
L e L 

<^) - Cos A1 

l 22 

Sin A1 SinAj* Sin (W g t^ + 

<f>^) T Cos A1 

23 

-Sin A1 Cos A^ 


31 

* . 

-Cos A1 SinA, Cos (W t. + 
L e L 

■ * L ) - Sin A1 

32 

-Cos A1 SinA, Sin (W^t, + 
L e L 

4>^) + Sin A1 

33 

* 

Cos A1 Cos A^ 



Where: 

★ 

;. L = Geodetic latitude of launch site 
<f> L = Longitude of launch site 
W o = Angular rate of the earth 
= Time of launch (from epoch) 

A^ = Launch azimuth 


Sin (W e 't L ..+ <J L ) 
Cos (W fi t L + 4> l ) 

Sin (W fi t L + d> L ) 
Cos (W e t L + f L ) 
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AVEH (Continued) 
ROTATION EQUATIONS 


2.1.4 
2. 1.4. 1.2 


DESCRIPTION 

This model defines the angular accelerations of the rigid body assuming 
that the center of mass lies in the X-Z plane (Iy z = 1^ = 0). 


MATH MODEL 
1. Equations 

q = *\ + pr (I zz “ ! xx^ + " p2 ^ ! XZ 

p = (a I zz + b ! XZ )/C 


r = (a I xz + b I xx )/C 

a = eM y + qr ^YY ! ZZ^ + pq ! XZ 
' X B 

b - eM z + pq (I xx ~ iyy) ~ q *^ ^XZ 
B 

2 

c = I T _ T c 

c XX L ll x XZ 


2. Definition of Symbols 

q = Angular acceleration about the Y body axis 

p = Angular acceleration about the X body axis 

r = Angular acceleration about the Z body axis 

*XX 5 *YY’ *ZZ = ^ oment °'* : i nert i a about X, Y, Z body axis respectively. 
! xz = x - z Cross product moments of inertia 
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2.1.4 AVEH (Continued) 

p, q, r = Integral of p, q, r (Body rates) 

zMy , zM v , eM 7 = Sum of moments about X, Y, Z body axes 
a B t B B 

= Aero moments + thrust moments + RCS moments 
+ engine deflection moments + slosh moments 


INPUT/OUTPUT 

Inputs: From THRUST, AERO, RCS, TVC and SLOSH 

Moments (about body axes) due to aerodynamics, main propulsion, 
reaction control, engine accelerations and slosh. 

Outputs - P, Q, R, P Q, R to IMU, Aero and Euler Angles 
IV. COORDINATE SYSTEMS 

Body axes - See translation equations for defi niton. 
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2.1.4 AVEH (Continued) 

2. 1.4. 1.3 EULER ANGLES 

I. This model defines the rate of change of the euler angles describing 
the attitude of the vehicle in inertial space. 

II. Math Model 

0 = (q cos <J> - r sin <>)/cos i> 

4 = qsin<{>+r cos 4 
$ = p = tan 4 (q cos 4. + r sin4>) 

Where: 

e, 4, 4> = Euler angle rates (1st, 2nd, and 3rd rotations, respectively) 
e, 4, 4 = Integral ofe, 4, 4 


III. Input/Output 

Inputs p, q, r from rotation equations 
Outputs 0, \p,<p to[B] and to IMU 

& 
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2 . 1.5 
2 . 1 . 5.1 
This model 
attitude i 
model . 


2 . 1 . 5.2 


DETERMINE 

COMMAND 

SIGNALS 


CALCULATE 

INERTIAL 

ATTITUDE 

ERRORS 


SCALE 
AND 
FILTER 
B SIGNAL 


INTEGRATE 

0 

SIGNAL 


BLCONT (Baseline Control System) 

PROGRAM DESCRIPTION 

issues commands to the engine gimbals such that the actual 
s made to follow the attitude required by the guidance 


MATH MODEL 
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2.1.5 


BLCONT (Continued) 


TRANSFORM 

TO 

BODY . 

ERRORS 

ADD 

0 INTEGRATOR 
TO ROLL ERROR 


SCALE AND 
FILTER 
ERROR 
SIGNALS 


FILTER 

AND 

BLEND 

ACCELEROMETER 

SIGNALS 


CALCULATE 

SCALE 

AND 

FILTER 

ACCELERATION 

ERRORS 


SUM 

FILTER 

AND 

LIMIT 

ERRORS 









2 . 1.5 


BLCONT (Continued) 


INTEGRATE 

ERROR 

SIGNALS 


TRANSFORM 
RATE 
COMMAND 
TO BODY 


FILTER 

AND 

BLEND 

RATE 

GYRO 

SIGNALS 


SCALE 

AND 

FILTER 

RATE 

SIGNALS 

ADD 

g SIGNAL TO 
ROLL RATE 
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5 


BLCONT (Continued) 


SUM 

ERROR 

AND 

RATE 

SIGNALS 


LIMIT 
ROLL AND 
RUDDER 
COMMANDS 


COMPUTE 

INDIVIDUAL 

ORBITER 

ACTUATOR 

COMMANDS 


LIMIT 

SRM 

SIGNALS 
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2.1.5 


BLCONT (Continued) 


INTEGRATE 

AND 

LIMIT 

SRM 

SIGNALS 


PREFERENCE 

ROLL 

COMMANDS 

OVER 

PITCH 

COMMANDS 


COMPUTER 

INDIVIDUAL 

SRM 

ACTUATOR 

COMMANDS. 
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2. 1.5.3 

SYMBOL DEFINITIONS 

ALIM 

Aero Error Software Limit 

KA (t) 

Aileron Roll Rate Gain 

KA' (t) 

Aileron Roll Error Gain 

KAY (t) 

Y Acceleration Gain 

KAZ (t) 

Z Acceleration Gain 

Kle 

Sideslip Feedback Integrator Constant 

KISRM 

SRM Integrator Constant 

KP (t) 

TVC Roll Rate Gain 

KP i (t) 

Individual Roll Rate Gyro Gain . 

KPA (t) 

Pitch Channel Error Gain 

KQ (t) 

Pitch Rate Gain 

KQ i (t) 

Individual Pitch Rate Gyro Gain 

KR (t) 

Yaw Rate Gain 

KR i (t) 

Individual Yaw Rate Gyro Gain 

KY i (t) 

Individual Y Accelerometer Gain 

KYA (t) 

Yaw Channel Error Gain 

KZ i (t) 

Individual Z Accelerometer Gain 

K0 

Pi tch Error Gai n 

K<j) 

Roll Error Gain 

K<|, 

Yaw Error Gai n 

PC 

SRM Pitch Error Signal 

PER 

Limited Pitch Channel Error 

PERR 

Integrated Pitch Channel Error 

PLIM 

Pitch Error Software Limit 
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2. 1.5.3 SYMBOL DEFINITIONS (Continued) 


PRATE'/ 

PRATEC 

PRATEF 

PRATEM 

PSRM 

RAERO 

RC 

RERR 

RLIM 

RRATE i 

R RATE A 

RRATEF 

RRATEM 

RSRM 

SRMLIM 

SRMRL 

V 

YC 

YF 

YLIM 

YM 

YER 

YERR 

YRATE i 

YRATEC 


Individual Pitch Rate Gyro Signal 
Pitch Rate Command 
Scaled Pitch Rate Error 
Blended Pitch Rate 
SRM Pitch Integrator 
Limited Roll Error for Aileron 
SRM Roll Error Signal 
Limited Roll Error for TVC 
Roll Error Software Limit 
Individual Roll Rate Gyro Signal 
Scaled Roll' -Rate for Ailerons 
Scaled Roll Rate for TVC 
Blended Roll Rate Signals 
SRM Roll Integrator 
SRM Integrator Limit 
SRM Error Signal Rate Limit 
Individual Y Accelerometer Signal 
SRM Yaw Error Error Signal 
Scaled Y Acceleration Error 
Yaw Error Software Limit 
Blended Y Accelerometer Signals 
Limited Yaw Channel Error 
Integrated Yaw Channel Error 
Individual Yaw Rate Gyro Signal 
Yaw Rate Command 


36 


2. 1.5. 3 SYMBOL DEFINITIONS (Continued) 


YRATEF 
Y RAT EM 
YSRM 

*1 

ZD 

ZF 

ZM 

e 

bint 

bt 


6 ip’ 6 iy 
<5D 


6PITCH 

6R0LL 

6RUDDER 


6SRM ip | 
6SRM iy f 
6YAW 


ec 

ec 

eEB 

eEi 

eF 


T 


Scaled Yaw Rate Error 

Blended Rate Gyro Signals 

SRM Yaw Integrator 

Individual Z Accelerometer Signal 

Z Acceleration Command 

Scaled Z Acceleration Error 

Blended Z Accelerometer Signals 

Vehicle Sideslip Angle 

Sideslip Feedback Integrator Signal 

Scaled and Filtered Sideslip Signal 

Individual Pitch and Yaw Engine Deflection Commands 

Prestored Engine Deflection Command 

Total Pitch Channel Engine Deflection Command 

Total Roll Channel Orbiter Engine Deflection Command 

Rudder Deflection Commands 

Individual Pitch and Yaw SRM Engine Deflection Commands 

Total Yaw Channel Orbiter Engine Deflection Command 

Platform Inner Gimbal Angle 

Platform Inner Gimbal Angle Command 

Platform Inner Gimbal Rate Command 

Body Pitch Attitude Error 

Platform Inner Gimbal Angle Error . 

Scaled Pitch Error 

Integrator Gain 

Platform Outer Gimbal Angle 
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2. 1.5. 3 

SYMBOL DEFINITIONS (Continued) 

<|>C 

Platform Outer Gimbal Rate Command 

<JrA 

Scaled Roll Error for Ailerons 

tpC 

Platform Outer Gimbal Angle Command 

4>EB 

Body Roll Attitude Error 

tj)EI 

Platform Outer Gimbal Angle Error 

<f>F 

Scaled Roll Error 

'P 

Platform Middle Gimbal Angle 

tpC 

Platform Middle Gimbal Angle Command 

\pEB 

Body Yaw Attitude Error 

ipEI 

Platform Middle Gimbal Angle Error 

'pF 

Scaled Yaw Error 


2. 1.5.4 INPUT/OUTPUT - 

Required inputs to this model are: 

• Inertial attitude angles 

• Body rotational rates 

• Y & Z translational accelerations 

• Prestored engine deflection commands 

• Prestored accelerations commands 

• Attitude commands from guidance system 
o Vehicle sideslip angle 

The system outputs are: 

o Deflection commands to orbiter engines 

• Deflection commands to SRM engines 

ft Deflection commands to aerosurfaces (rudder) 
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2.1.6 CGAINS (Control Gains Equations) 

2.1.6. 1 PROGRAM DESCRIPTION 

The CGAINS program is used to calculate the control gains necessary 
for a desired type of control during the Shuttle boost. There are 
several options for the control gains that are calculated: load 

minimum, drift minimum, or attitude control for the pitch and yaw 
gains; and thrust vector control or aileron control for the roll 
gains. The following model presents the equations necessary to 
calculate these gains. 


2. 1.6. 2 MATH MODEL 


The following quantities must be calculated each time the control 
gains are needed. The symbols used in these equations are defined 
in Table 1. 


Ip 


iy 



<W C zO> + X cg 


(C /C ) + X 
' ma' zor eg 


< b W + * 


’ X R 
X R 

eg “ X R 


N 'p 

= qS C 

Za 

% 

= q S C 

H ye 

N o 

= q s c 20 

] ap 

= -i 

ap 

1 

ay 

= -1 

ap 

f ax 

= q S C x0 
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1.6 


CGAINS (Continued) 


'Ip 


■ (F 0 cos 6 CO + F B cos S CB + r k )M 


K 2p N'p/m 


K 


3p 


= -Fg COS 6gg/m 


K. = -K,_ 
iy ip 


K 


2y 


= -N' y /m 


K 3y = 2/3K 


3p 


= Ng/m 


■' X cg - X o 


V = Z cg' Z 0 


lg =(’ x 2 + l z 2 )" 


e 0 ( ] o -/AX ^z^p 


6 = tan 

CGO 


-1 -V 


C 2p Vo C0S Vo " ^CGO^p 


A = C„ K„ - C 1 IC 
p 2p 2p lp 3p 
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'2 y 


' ,; ■ - - •• ■• **>/ Va-»- “* .-, i- j 

= (2/31 x) F/I , 


.• " • ' ■ • * ,♦ 
;* .' ^ V‘* ■ 'r\s 
■ * - ^ ^ . ’’ * 

•.*: . • 

_• .■ 

■ . - '; . 

' v 

v - v.?* ■•* 


. ? = • j&. \ •’■ , , . 




. : c 2y K 2y ; , C ly K 3y ; ; 


•Vi/' 


2.1.6. 2.1 PITCH AND YAW CONTROL GAINS FOR LOAD 




MINIMUM OPTION 


x-aMA-WtS' _ : „ 


a Op = °- 


3 2p 


0). 


- C 


ilL 


A + cu 

p y 


(K 3p + % C 2p) 


IP 


iipJjL 

C 2p 


1 g 2p ^ K 3p + ^ap C 2p^ 


a 0y = °* 


2y 


«, - C. 

-i Ll 


A„ + », (K, + 1 Co ) 
y z v 3y 'ay 2y ; 


a, = 


ly 


'2y 


1 “ g 2y (K 3y + ] ay C 2y } 


41 


2.1.6 CGAINS .... (Continued) L 

2.1. 6. 2.2 PITCH AND YAW CONTROL GAINS FOR DRIFT MINIMUM OPTION 



“v - * C 2p WV C lp 


Co + A +03 

2p lp p y 


< K 3 P + V V 


a 0p g 2p K lp + C lp K lp /A p 


ip c 2r 


2 c to 
T_.il 


1 - g 2p + >ap C 2p } 


2y 


“ Z " ^ + c 2y K ly /r \^ C ly 


C 2y K ly + A y + “z (K 3y + ] ay C 2y ) 


a n = go K. + C, K. /a 
Oy y 2y ly ly ly' y 


a. - 


iy c 


2y 


1 - g 2y < K 3y + V C 2y> 
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2.1.6 


CGAINS 




2. 1.6. 2. 3 PITCH AND YAW CONTROL GAINS FOR AniTUDE CONTROL 


• _4 ■ ■ r r? ■ ■= 

; >'< v >-S’ - 


g 2p 


; = 0 


(o C-i 

a o P = y lp 

C 2p - 


‘IP 


2 c to 

■ p y 

C 2p 


g 2 y 


= 0 


A 


x 0y 


w - C-. 

= z ly 


"’*■ C 


2 y 


‘iy 


2 C u) 

y z 


'2y 
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2.1.6 CGAINS (Continued) 

2. 1.6. 2. 4 ROLL CONTROL GAIN EQUATIONS 
2. 1.6. 2. 4.1 Thrust Vector Control 


c 2r - VC + Z o 2 ' F/I r 

2 

OJ 

a Or = **-■ 


'2r 


lr 


2 c r 
C 2r 


2.1. 6. 2. 4. 2 Aileron Control 


q S b C-j 


'2r 


w x 2 


Or C 


2r 


2 c u) 

x 


lr C 


2r 


2. 1.6. 2. 4. 3 TVC and Aileron Control 


'2r 


Vv 2 + Z 2 F + q S b C, 
» o o . 1 


ag r and a^ r same as above 
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2 . 1,6 _ 
2 . 1 . 6. 3 
VARIABLE 


3 1P 

a iy 


c 

c 

C 

C 

C 


mo* 

zo’ C za 


ye 

C nB 

% 

C xo 


1 B 
F AX 
g 2 P 

9 2y 

I 

K 

1 

m 


N' 


q 

S 


CGAINS (Continued) 

DEFINITION OF SYMBOLS 

DEFINITION 

Pitch attitude gain 

Yaw attitude gain 
Pitch attitude rate gain 
Yaw attitude rate gain 
Wing span 

Coefficients for calculation of gains 

Mean aerodynamic chord 

Aero, moment coefficients 

Aero, normal force coefficients (pitch) 

Aero, normal force coefficients (yaw) 

Aero, moment coeff. (yaw) 

Coeff. of roll moment with aileron 
Aero, axial force coeff. 

Total orbiter thrust 
Total booster thrust 
Aero, axial force 
Pitch acceleration gain 
Yaw acceleration gain 
Moments of inertia 

Coefficients for calculation of gains 

Moment arms 

Mass 

Normal force at zero angle-of-attack 
Partial of normal force 
Dynamic pressure 
Aerodynamic reference area 
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2.1.6 CGAINS (Continued) 

VARIABLE DEFINITION 

X Q , Y , Z Q Average location of orbiter engines 
X CG* Y CG Location of center of gravity 

1 Q Translational acceleration at zero angle-of-attack 



Temporary variable 
Damping ratio 

Rotational acceleration at zero angle-of-attack 
Natural frequency 


2. 1.6.4 PROGRAM FORMAT, INPUT-OUTPUT 

The control gains routine calculates the attitude gains, attitude 
rate gains and accelerometer gains for each of the desired condi- 
tions which have been previously mentioned. These gains are also 
calculated for a range of natural frequencies (co). These results 
are output on a scratch tape or FASTRAND file for processing by a 
plotting program. 
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2.1.7 


ACTIVE GUIDANCE 


2. 1.7.1 PROGRAM DESCRIPTION 

Active guidance is used to provide inertial steering commands during 
the boost to orbit insertion. Described in this section are two active 
guidance schemes which are incorporated in the SSFS. These two guidance 
schemes are the Apollo Lunar Ascent Guidance or 1 E ' Guidance (program 
name - BACTGD);and Linear Tangent Guidance (program name - GUIDE). 

2. 1.7.2 BACTGD 

The 'E 1 guidance method is baselined in Revision B. of. "MSC Space Shuttle 
GN&C Design Equations Document". The math model for this method is described 
as follows : 


Preliminary Calculations : . ' 

R p (4) = (Rp(l) 2 + R p (2) 2 + R p (3) 2 )^ 

G ( 1 , 1 ) = Rp(l)/R p (4) 

G(1 >2) = R p (2)/Rp(4) 

G(l,3) = Rp(3)/R p (4) 

G ( 3 , 1 ) = 6 ( 1 , 2 ) U Q ( 3 ) - G(1 ,3) U Q (2) 

G (3 ,2) = G(1 ,3) U Q (1) - G( 1 ,1) U Q (3) 
G(3,3) = G(l,l) U Q (2 ) . - 6(1,2) U Q (1) 

G ( 2 , 1 ) = G(3 ,2) G(1 ,3) - G(3,3) G(l,2) 

G(2,2) = G(3,3) G(l,l) - G(3,l) G(l,3) 

G(2 ,3) = G(3,l) G(1 ,2) - G(3,2) G(l,l) 
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2.1.7 GUIDANCE (Continued) 


VU' 


v p (if 

V 2 > 

- [e] 

V p (2) 

_V S (3)_ 


V p (3) 


v gd (D = R - v G d) 

V GD (2) = V - V g (2) 

V GD ( 3 ) « Z - V g (3) 

V PD t G ] V GD 

AM(1) = G(l,2) V p (3) - G(1 ,3) V p (2) 

AM(2) = G ( 1 , 3) V p ( 1 ) - G(l,l) V p (3) 

AM (3) = G(l,l) Vp ( 2 ) - G(l,2) V p ( 1 ) 

AM(4) = Vam(1) 2 + AM(2) 2 + AM(3) 2 * 

a grav^ -vW>- + a grav^ + a grav^ 

G EFF = " A GRAV (4) + AM ( 4 )/ R p( 4 ) 

V PD (1) = V pD (l) - %T G0 G EF|r G(l,l) 

v pd {2) = W 2) - ^ t go g eff G(1 ’ 2) 

v p D (3) = V po^ 3 ) " ^ G0 G^pp G(l,3) 

V PQ^ 4 ) ~ q (^) ~ VpQ (2) + v po^ 3 ) 


NOTE: Prior to the end of Boost Polynomial Guidance the "Preliminary 

Calculations" portion of Active Guidance shall be cycled 5 
times to establish an initial value of T G q using the following 
equation: V p[) (4) 
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2.1.7 GUIDANCE (Continued) 

CALCULATIONS FOR CONSTANT THRUST: 
V r 


T TL t " A 


EG 


TL 


Qj = In (1 - -^) 


i 


t go t tl + 


v p D ( 4 ) + V EG Q 1 

a tl 


Q 2 “ V EG Q 1 
Q 3 = T GO " T TL 

A 11 = <z + ^3 a tl 

Q 4 = tQ 2 
Q 5 = V EG T TL 

A 12 = " Q 4 ~ Q 5 + ^TL ^ T GO ? - T 


TL 


A 21 Q 4 + Q 5 " Q 2 T G0 + JsA TL Q 3 


a 22 - T - <«4 + Os - We) + 5 <-r - T ' 

+ — (T ^ - 3T T ^ + 2T 
6 GO’ J GO *TL TL ; 


GO 


CALCULATIONS FOR CONSTANT ACCELERATION: 


T 


V po (4) 


GO A tl 
A = A T 

M 11 h tl ‘go 

A 12 = ^11 T GO 


A 21 A 12 

A 22 = 1/3 A 12 T GO 
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2.1.7 GUIDANCE (Continued) 

FINAL CALCULATIONS: 

V = A n A 22 - A 12 A 21 


Q 2 = U Q (1) R p ( 1 ) + U q { 2) R p (2) + U q ( 3) R p (3) 


A = 


-A 12 (R - R p (4) - R T go ) + A 22 V g (1) 


B = 


A n (R - R p (4) - R T go ) - A 21 V G (1) 



r _ ' A 12< V - 02 - V W +A 22 V 2 > 

V 

D _ A 11< Y - g 2 - j W - A 21 V 2 > 

..... .. V 

A g (1) = I A T I (A + BcSx) - G eff 
A g ( 2) = |A T |(C + D 6 x ) 

% = |a t | 2 - a g (i) 2 - a g (2) 2 

A g (3) = AMAX1(0 . , VO^) 


Q 7 = a g (3) 2 + A g (1) 2 


'V -■ W T a g - 

«p(4) = Y A p a) 2 + Ap(2) 2 + Ap(3) 2 
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SSFS MODEL DOCUMENTATION SERIES 


2.1.7 GUIDANCE (Continued) 



i 

i 

Symbol Definitions’ 




h 

U Q 


V G 

0 0 «« 
R,Y,Z 

V GD 

VpD 

AM 


GRAV 

*EFF 

f G0 

; EG 


A 1 

T 

T 


TL 


/ 

transformation from platform coordinates to guidance 
coordinates 

current position vector in platform coordinates 
target unit vector normal to orbit plane in platform 
coordinates 

current velocity vector in platform coordinates 
current velocity vector in guidance coordinates 
target, veloci ty in guidance coordinates 
velocity to be gained in guidance coordinates 
velocity to be gained in platform coordinates 
angular momentum- vector 

acceleration due to gravity in platform coordinates 

effective gravity 

time to go until orbit insertion 

rocket exhaust gas velocity 

acceleration due to thrust in platform coordinate = 

M 

time required to burn up the vehicle's total mass 
time until thrust limiting 
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2.1.7 


GUIDANCE (Continued) 


Ay L = vehicle acceleration at the acceleration limit 

Qj. Q 2 , Q 3 .Q 4 .Q 5 , Q 6 , Q 7 = temporary variables 

XCG, YCG, ZCG = coordinates of vehicle center of mass 
ZCF = position of thrust centerline 

A,,, A 19 , A 91 , A„„ = parameters needed to compute linear control 
11 1 1 ^ coefficients 

V = temporaty variable 

A, B. C. D = 1 linear control coefficients 

6t = time period between guidance calculations 

Ag = desired acceleration in guidance coordinates 

7tp = desired acceleration in platform coordinates 

e D1 = temporary variable 

e D’ ^D’ ~ desired platform gimbal angles 

R, Y = target position in guidance coordinates 

2. 1.7. 3 GUIDE 

The LTG guidance method is baselined in Revision C of "MSC Space Shuttle 
GN&C Design Equations Document" , and also in NASA Internal Note MSC-IN-72-39. 
The LTG math model is as follows: 


GUIDANCE (Continued) 



EXIT TO GUIDANCE 


DESIRED ORBIT 

COORDINATE 

SYSTEM, 

G MATRIX 


A^ = Launch Azimuth 

= Geodetic Launch Latitude 

i Q = Desired orbit inclination 

0p = Desired Longitude of descending 
node (measured from launch 
meridian) 





2.1.7 GUIDANCE (Continued) 



INITIAL RANGE ANGLE 


LOCAL GUIDANCE 
COORDINATE SYSTEM 


RADIUS AND VELOCITY IN 
LOCAL GUIDANCE FRAME 

MEASURED ACCELERATION 
MAGNITUDE 


AVERAGE RADIAL GRAVITY 
MAGNITUDE 
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GUIDANCE (Continued) 


ENGINES 

ON 


k= NO. OF PRESENT BURNING STAGE 


T BK -■ 0 

t ck ■ V T 



T k = V exk /a 


T C ■ T+T Ck 


?T ?o 



FIRST 

PASS 


RANGE ANGLE-TO-GO IN 
DESIRED ORBIT PLANE 


$GO * 5 [( V G3/ R + ( v j) cos Y D )/ R D ] T gq 


V T1 ’ - V D sin *4»S0 - V 


V T2'° 

TERMINAL 

V T3 * V D C0S (< * ) GO - V 


R 71 = R D cos< f , GO 

VELOCITY AND RADIUS VECTORS 

R T2 " 0 

IN LOCAL GUIDANCE FRAME 

R T3 = R D sln( ^Ga 



V 







2. 1.7 GU I DANCE ’ " (Conti nued ) 



GRAVITY MODEL 


AVERAGE GRAVITY VECTOR 
FOR VELOCITY LOSSES 

AVERAGE GRAVITY VECTOR . . 
FOR DISTANCE LOSSES 


ACCELERATION LIMIT 
LIMIT TO BURN OUT 
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2.1.7 GUIDANCE ' (Continued) • • 





2.1.7 GUIDANCE (Continued) 



TIMES USED FOR MULTI-STAGE 
INTEGRALS 


K = No. of present burning stage 
n = No. of final stage 


T RUi 

T 0Ui 


T RLi * T LSi 


E 

j = k 


T Bj - T Bi 


T 0Ui + T Li 



Time remaining after constant- 
acceleration Durn 

Time remaining after constant 

thrust burn 

Initial Time of constant thrust 
burn 

Initial Time of constant 
acceleration burn 


i = i + 1 {4 


L = 
J = 
S = 

Q = 

0 

0 

0 

0 

J 

N 

L 

i = 

DU 


Initialize Summation Loop 


V7 
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2.1.7 GUIDANCE (Continued) 
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, _ 2 .1; 7- - GUIDANCE - (.Continued) - 



OISTANCE-TO-GO 



(Down range component) 



ACCURATE RANGE ANGIE 
FOR NEXT GUIDANCE PASS 









2.1.7 GUIDANCE (Continued). 



EXIT MAJOR CYCLE LOOP 


STEERING VECTOR RATE 

TRANSFORM COMPONENTS OF 
STEERING VECTOR FROM 
LOCAL GUIDANCE FRAME TO 
PLATFORM FRAME 
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2.1.7 


GUIDANCE (Continued) 

2. 1.7.4 INPUT-OUTPUT 

e e * 

The targeting program provides lL, R, Y, R, Y, and Z. The EOM 
program provides Rp, Vp and A^. The gravity program provides Ag R ^y 
The resulting quantities calculated by the model aree^, and <j,p. 

No flight software commands are accepted by this model. The EOM, 
Targeting, and Gravity programs must be present to provide inputs to 
the model . 

2. 1.7. 5 COORDINATE SYSTEMS 

The three systems used by this model are the Guidance, Platform, and 
Body systems. The center of the guidance system is at the earth center 
of mass. The R axis points from the center of the earth to the vehicle, 
the Y axis is perpendicular to the orbit plan pointing to the pilot's 
right, and the Z axis completes the right-handed triad. The body system 
is fixed with respect to the vehicle with the X axis forward through the 
main propellant tank centerline, the Z axis in the engine. gimbal pivot 
plane and pointing down, and the Y axis points toward the pilot's right 
completing the right-handed triad. The platform system origin is at 
the earth's center of mass and is fixed in inertial space at the time of 
launch. The X axis is parallel to but opposite in sense to the launch 
pad gravity vector, the Z axis points downrange in the launch plane, and 
the Y axis points toward the pilot's right completing the right-handed 
triad. 
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2.1.7 GUIDANCE (Continued) 


to- 


w ” 


i*v- 


cose 

0 

sine 

0 

1 

0 

-sine 

0 

cose 

! 

j 

j 

CO Sip 

- 

i 

si nip 0 

si nip 


cosip 0 

0 


0 1 

I 

1 

0 

i 

j 

- — 

0 

0 

COSip 

-sin<p 

0 

sin<p 

COSip 


Xp 


'v 1 
X B 

V P 

• i-i [*] i=; 

y b 

Z P. 


h. 


where Xp, Y p> Zp = platform coordinate vector 
Xg , Yg, Z B = body coordinate vector 
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2.1.8 HINGE (HINGE MOMENT CALCULATION) 

\ 

2. 1.8.1 PROGRAM DESCRIPTION 

This program returns the rudder hinge moment coefficient as a function 
of rudder deflection, angle of attack, sideslip angle, and Mach number 
This is achieved by a table lookup for hinge moment coefficients due 
to sideslip ( Ch ^ ) and hinge moment coefficient due to rudder de- 
flection angle (Ch f ), each as a function of angle of attack and Mach 

R 

number. Tables of CH n (M) and Ch. (M) for a = 0°, 5°, and 10° are 

6 S R 

used. If the angle of attack is greater than 10°, 10° is used for - 
lookup; if angle of attack is less than 0°, 0° is used. A double 
linear interpolation is then used for CH„ and Chi- -as a function of 
angle of attack and Mach number. 


2. 1.8.2 MATH MODEL 

DATA / Ch^ (M)a = 0 / 

DATA / Ch g (M)a = 5 / 

DATA / Ch g (M)a = 10 / 

DATA / Ch 6R (M)a = 0 / 

DATA / Ch fiR (M)o = 5 / 

DATA / Ch 6R (M)a = 10 / 

IF a >10° a =10° 

IF a <0° a = 0 ° 

TABLE L00PUP AND INTERPOLATION FOR, 
Ch g = Chg (M, a) 

Ch. _ Ch (M, a) 

Ch = Ch + Ch * 6 d 

P 0 R K 


64 


2 . 1.8 


HINGE (HINGE MOMENT CALCULATION) Continued 


2. 1.8.3 INPUT/OUTPUT 

Formal parameters in the call to HINGE are 1) rudder deflection 
angle (degrees), 2) angle of attack (degrees), 3) sideslip angl 
(degrees), 4) Mach number, 5) hinge moment coefficient to be 
returned, and 6) KDATA, a provision for using the routine for 
different aero surface coefficients, KDATA = 1 for rudder. 



2.1.9 MASPRO (Mass Properties) 

2. 1.9.1 PROGRAM DESCRIPTION 


This model provides the mass properties, which consist of center of 
gravity travel, moments of inertia as a function of weight and total 
mass calculation., 

2. 1.9.2 MATH MODEL 


w 

= 

table lookup 

F(H, 

M 

= 

M - W AT/gc 


IF Ml 

> 

M > M2, M = M 

- (Ml 

W 

= 

Mg c 


CG x 

= 

table lookup 

F(W) 

CG z. 

= 

table lookup 

F(W) 

*xx 

= 

table lookup 

F(W) 

l yy 

=" 

table lookup 

F(W) 

hz 

= 

table lookup 

FfW) 

! xz 

= 

table lookup 

F(W) 


Where: 

W = total vehi cle weight 

M = total vehicle mass 

Ml = mass prior to abort SRM jettison 

M2 = mass after abort SRM jettison 

W = weight flow rate 

g c = mass to weight conversion constant 
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2.1.9 MASPRO (Mass Properties) Continued 
aT = computational time increment 

CG = X center of gravity location from vehicle reference 

A 

CG Z = Z center of gravity location from vehicle reference 

I = Mass moment of inertia about X axis 

XX 

I = Mass moment of inertia about Y axis 

I = X - Z cross product of inertia 

H = Altitude 

P = Fractional throttle setting 

2. 1.9. 3 INPUT/OUTPUT 

Input g c and initialize M and H from cards 
Read W, CG and Moments of Inertia from tables. 

Output M, CG's and I 's. 
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2.1.10 


MAXMIN (Maximum and Minimum Values Printout) 

2.1.10.1 PROGRAM DESCRIPTION 

The MAXMIN program is used to printout the maximum and minimum values 
of several flight parameters at any time that it is scheduled. At 
any time that MAXMIN is called, the maximum-minimum values accumulated 
until that time is printed. The program may be called as often as 
desired. 

2.1.10.2 MATH MODEL 

;ln MAXMIN, all that is done is a series of WRITE statements. The 
actual saving of the desired maximum-minimum values and the times they 
occur is done in BLCONT. This is done in the form: 

it x > ia| , j * ; f F k > 

where X is the present value of the variable in question, A is the 
maximum value of the variable, T is the time that the maximum occurred, 
and t is present time. Although this is done in BLCONT, these values 
could be saved in any routine that is called every program cycle. 

2.1.10.3 INPUT/OUTPUT 

All of the variables which are output by MAXMIN are transferred 
from BLCONT in common blocks A^ x and T^. The following quantities 
are output: 
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2.1.10 MAXMIN (Maximum and Minimum Values Printout) Continued 


q 

= 

dynamic pressure 

a 

= 

angle-of-attack 

6 

= 

sideslip -angle 

q-a 


Q-Alpha 

q-e 

= 

Q-Beta 

p 

= 

roll (inertial) 

p 

= ' 

roll rate 

6 a 

= 

ai leron deflection 

*a 


- aileron deflection rate 

V 

= 

rudder deflection 

K 

= 

rudder deflection rate 

6 e 

= 

engi ne deflections 

r hm 

= 

rudder hinge moment 

a hm 

= 

aileron hinge moment 

r hp 

= 

rudder horsepower hours 


= Foonman's Constant (equals 1.0) 
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2.1.11 ORBITI (Orbit Parameters) 

2.1.11.1 PROGRAM DESCRIPTION 

This math model calculates the parameters of the trajectory achieved 

" . ; 

at insertion from knowledge of the state vector in polar equatorial 
coordinates at the time of orbiter engine shutdown. The program 
calculates node, inclination angle, orbit phase angle, eccentricity, 
orbit parameter, true anomaly, apogee altitude, and perigee altitude 

2.1.11.2 MATH MODEL 


n = Arctan 


~ V YPE R ZPE + V ZPE' R YPE 
-V XPE R ZPE + V ZPE R XPE 


C = Arctan 


_ZPE 

" V XPE Sinn + V YPE ' C0Sn 


) 


v = Arctan 


-Vxpe cosn — Vyp E si nrj 
cos ? (" v xp£ sinn + v yPE cosn ') + V zpE sin 




2.1.11 ORBITI (Orbit Parameters) Continued 



TA = Arctan ^ 


R p = SMA ( 1-E) - R r 
R a = 2 • SMA - R p - 2 R e 


2.1.11.3 SYMBOL DEFINITIONS 

V XPE’ V YPE’ V ZPE = Velocl ' 1: Y in polar-inerti al 

R XPE’ R YPE’ R ZPE = Position in polar-inertial 

n = Longitude of the ascending node 

C = Inclination angle 


coordinates 

coordinates 
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2.1.11 ORBITI (Orbit Parameters) Continued 
e = Eccentricity 

K = Gravitational constant = 1.407654 x 10^ ft'Vsec^ 
v = Orbit phase angle (angle between equator and perigee) 

P = Orbit parameter 

SMA = Semi-major axis 

FP = Tangent of flight path angle 

TA = True anomaly 

Rp = Altitude at perigee 

R q = Altitude at apogee 

R e = Mean radius of the earth 

2.1.11.4 INPUT/OUTPUT 

The ORBITI math model requires the vehicle state vector in platform 
equatorial coordinates as input. As output the model prints out 
ascending node, inclination angle, orbit phase angle, eccentricity, 
orbit parameter, true anomaly, apogee altitude, and perigee altitude 
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2.1.12 ORB I fR (3D EQUATIONS OF MOTION) 

2.1.12.1 PROGRAM DESCRIPTION 


This Math Model calculates vehicle acceleration from active guidance 
commands and integrates to get velocity and position in platform 
coordinates. Polar-Equatorial and Local Vertical coordinates 
systems are erected to calculate latitude, longitude, and flight path 
angle. Logic is included for acceleration limiting, integration 
cycle time rectification, and velocity cutoff. 

2.1.12.2 MATH MODEL 



N 

I 

1 = 1 


T v 0) 



<AP 


A F 

1 *; 


A 1 

m pt 


- Q pt p t Q ap + a grav 










2 . 1.12 


ORBITR (Continued) * '■- : 

“V l -' ’ :r ' s 5;'-j£'f--\P'£'2&£ 


V 


- A A 


PT dt + v r 


R^ 1 

R P 


■ /t 


t V p dt + R p 


W T R p ' 


= Arcs in 


Rp{3) 


= Arctan 


R F (2)' 

R^ID 


“e ^L +t > 


01 


V LV 


T W T V'p 


. A, V < 2 M ’ 

Arcs in / -jw — 

'LV 1 


Integration Cycle Time Rectification 

Acceleration limiting will most likely occur 
points. Therefore the integration for this 
two parts. , 

Velocity Cutoff 

The program will be terminated either when |’ 
is less than M vq , whichever occurs first. ~ 
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between integration time 
interval must be done in 


* 

p| exceeds Z or when My 



2.1.12 ORBITR (Continued) 
2.1.12.3 SYMBOL DEFINITIONS 


Qp T = Thrust acceleration at full throttle 

N = Number of engines 

T v = Maximum vacuum thrust of main engine 

M'y = Current vehicle mass 

\lM = Maximum allowed vehicle acceleration 

|A t | = Magnitude of acceleration due to thrust 

P T = Throttle setting 

My = Total mass flow rate for all engines at full throttle 

Ap = Acceleration command from guidance in platform coordinates 

Q Ap = Unit vector acceleration command 

A pT ' = Current total vehicle acceleration in platform coordinates 
Aqrav = Acceleration due to gravity in platform coordinates 

Vp = Velocity of vehicle in platform coordinates on last pass 

Vp' = Current vehicle velocity 

Rp = Position of vehicle in platform coordinates on last pass 

Rp' = Current vehicle position 

*-> 

Rp = Position of vehicle in polar equatorial coordinates 

•k 

Ap = Geodetic latitude of launch site 

* 

4>P = Longitude of launch site 

w = Rotation rate of earth 

e 

Tp = Time of launch 

t = Elasped time since launch 

= Mass of vehicle on last pass 


2.1.12 ORBITR (Continued) 

A^ = Launch azimuth 

Ay = Present vehicle latitude 

<j>‘ = Temporary variable 

[а] = Transformation from polar equatorial to platform coordinates 

[ б ] = Transformation from local vertical to polar equatorial 

coordinates 

<t> = Present vehicle longitude 

y = Flight path angle 

V^y = Velocity in local vertical coordinates 

Z = Target velocity in guidance coordinates 

My q = Mass of empty orbiter 

2.1.12.4 INPUT/OUTPUT 

The targeting program provides Z. The Active Guidance program 

provides A p . The gravity program provides A The resulting quan- 

titles calculated by this model are Ay, 4 ., y, V p , R p , A pT> My, and V p y. 

• ^ 

Flight software commands accepted by this model are A p . The Active 

Guidance, Targeting, and Gravity programs must be present to provide 

inputs to the model. 
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2.1.12 


ORBITR (Continued) 



The three systems used by this model are the Platform, Polar Equatorial, 
and Local Vertical systems. The platform system origin is at the 
earth's center of mass and is fixed in inertial space at the time of launch. 
The X axis is parallel to but opposite in sense to the launch pad gravity 
vector, the Z axis points downrange in the launch plane, and the Y axis 
points toward the pilot's right completing the right-handed triad. The 
polar equatorial system origin is at the earth's center of mass and is 
fixed in inertial space at the time of launch. The X axis is in the 
equatorial plane pointing toward a reference meridian, the Z axis points 
through the north pole, and the Y axis completes the right handed triad. 

The local vertical system origin is at the earth center of mass. The X 
axis .points from the earth center to the vehicle, the Z axis is in the 
plane containing the earth's rotation axis and the X^ v axis. The Z axis 
is perperdicular to the X axis and points toward the north pole. The 
Y axis completes the right handed triad. - 


[a] 



Rp = [a] Rp 


r* 

a ll 

a 12 

a is1 

a 21 

3 22 

a 23 

a 31 

a 32 

3 33 


77 



2.1.12 ORBITR (Continued) 


a ll 

= 

cos 

* 

X L COS U L 

* + Vl> 









* 

* 






a 12 

~~ 

cos 

Aj. sin (_* L 

+ u> e t L ) 






a 13 

= 

sin 

★ 










★ 


*. 





a 21 


sin 

A l sin Ap 

COS (u e t L 

+ 4> l ) 

- COS 

Ap sin 

(» e t L 

+ *p) 




* 


*. 




*. 

a 22 

= 

sin 

A l sin a^ 

sin (a» e t L 

+ 4*p ) 

- COS 

Ap cos 

(“e*! 

+ 4>p ) 

a 23 

= 

-si n 

★ 

1 Ap cos Ap 










★ 






* 

a 31 

= 

-cos 

Ap sin Ap 

COS (to t. 
e L 

+ 4>p ) 

- sin 

Ap sin 

(“ e t L 

+ ^ 




★ 


*. 




* 

a 32 


-cos 

Ap sin A^ 

sin (w t, 
e L 

+ 4>|_ ) 

+ sin 

A^ cos 

( "‘e\ 

+ 

a 33 

= 

cos 

-k 

Ap cos Ap 








-+ 

Rp = [s] Rj_y 


11 

d 12 

d 13 

21 

d 22 

d 23 

31 

d 32 

d 33 


sin 

A v 

= R f ( 3)/ |R f | 
r. 9 ' 

cos 

A v 

= Ml - sin^ a v 

si n 


= R F (2)/( | Rp 1 cos Ay) 

cos 


= Rp ( 1 )/ ( I Rp I cos A v ) 

d n 


= cos Ay cos <j>‘ 

d 12 


= -s i n <j) ’ 

d 13 


= -sin Ay cos <}>' 

d 21 


= cos Ay sin (j)' 
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2.1.12 ORBITR (Continued) 


= cos <j> ' 

d^^ = -sin Ay sin <j>' 

d^j = sin Ay 



d^2 COS Ay 
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2.1.13 


2.1.13.1 


ORBTAR (Boost Orbit Insertion Targeting Model) 
PROGRAM DESCRIPTION 


The Targeting program is used in the flight software to describe the 
orbit plane with respect to the launch pad. Position, velocity, and a 
unit vector normal to the orbit plane at perigee are. calculated from a 
knowledge of perigee and apogee altitudes, location of the launch pad, 
orbit inclination, and an orbit parameter. 


2.1.13.2 MATH MODEL 


Rn = 


R e + h P 


R A = 
A = 


R e + h A 

R p + R a 


R 

Y 
Z 

R 

Y 
Z 



= 0 

= unconstrained 


= 0 
= 0 
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2.1.13 


ORBTAR (Continued) 


S G = sin (x) cos(3)+ cos (A z ) cos (x) sin (3) 

C L = sin (A z ) cos(x) 



> = cos (y)/Cq 

A = Arcsin (Cp) -Arcsin (a) 

Ug (1) = -sin (e) sin (a) 

Up (2) = cos (A) 

Up (3) = cos (6) sin (A) 


Where : 

Rp = distance from earth center of mass to periapsis 

R^ = distance from earth center of mass to apoapsis 

h p = altitude at periapsis 

h A = altitude at apoapsis 

A = semi -major axis 

e = eccentricity 

R g = radius of earth 

R = radial distance at insertion 

Y = cross-range distance at insertion 

Z = downrange distance at insertion 

R = radial rate at insertion 
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2.1.13 


ORBTAR (Continued) 

Y = lateral velocity at insertion 

Z = downrange velocity at insertion 
Sg, C^, Cg, Cp, a, A = temporary variables 

X = latitude of launch pad 
u = universal gravitational constant 
3 = orbit parameter 

A z = launch azimuth 
y = orbit inclination angle 

— V 

IL = unit vector norma) to desired orbit plane in platform 
^ coordinates (see Section IV) 

2.1.13.3 INPUT/OUTPUT 

This model requires R g , hp, h flJ A , 3, A^, and Y as input and calculates 
R, Y, Z, R, Y, 1 , and Up. The model needs to be called only once per 
simulation. 

2.1.13.4 COORDINATE SYSTEMS 

This model calculates a unit vector in platform coordinates. The platform 
system origin is at the earth's center of mass and is fixed in inertial space 
at the time of launch. The X axis is parallel but opposite in sense to the 
launch pad gravity vector, the Z-axis points downrange in the launch plane and 
the Y axis points toward the pilot's right completing the right-handed triad. 
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2.1.14 RCS (Reaction Control System) 

2.1.14.1 PROGRAM DESCRIPTION 

This program calculates the moments and linear accelerations applied 
to the vehicle due to RCS thrust commands from the flight software. 

2.1.14.2 MATH MODEL 

Thrust commands are conditioned according to the following limitations 

(1) A jet cannot be commanded to ignite unless the duration of 
ignition is some minimum value 

(2) No jet can be ignited continuously longer than Tg seconds 

(3) A jet cannot be commanded to ignite unless T^ seconds has 
elapsed since the previous ignition has ceased. 

Figure 2-1 shows a logic diagram which can be used to implement these 
limi tations. 
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RCS LOGIC DIAGRAM 
FIGURE 2-1 
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2.1.14 RCS (Continued) 
i = engine index 

j = 1, 2, or 3 for x, y, or z respectively 
T MAX = ^rust achieved with continuous ignition 

T MIN = avera 9 e thrust for a computation cycle for minimum thrust duration 
TC^(j) = commanded thrust 
T;(j) = realized thrust 

= indicator for continuous thrust for the previous pass' 

(set equal to zero for restart) 

IN^ (j) = thrust enable flag j ■ . 

(set equal to zero for restart) 

T ch = check time for 

= check time for Tg 

t = current time at entry to RCS program 

n. = pi X 1\ (1) 

P.(l) = ELX. - XCG 

P i (2) = ELY i - YCG (2) 

Pi (3) = ELZi = ZCG 

^ . L. 

Pi = engine position vector for the i 1 jet; ELXi, ELYi, ELZi 

X, Y, and Z locations, respectively, of the i^ RCS jet cluster; 

XCG, YCG, ZCG = X, Y, and Z locations, respectively of the vehicle 
center of mass 
i = engine index 

— th 
Mi = moment vector due to thrust from the i jet 
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2.1.14 


RCS (Continued) 


T\ = thrust vector for the i t * 1 jet 
X = indicates vector cross product 


N 

— v T -*■ 

MR = • V i M i 


( 3 ) 


MR = total effective moment from all RCS jets 
N = number of RCS jets 

N 

F = .S j 
r i = 1 i 

F = total effective linear force 


( 4 ) 


A = reduction in mass of vehicle due to RCS fuel usage 
R = RCS jet efficiency constant 

2.1.14.3 INPUT/OUTPUT 

The model requires ELX, ELY, ELZ, XCG, YCG, ZCG, T , t, and M as input. 

c 

The flight software command accepted is T c< The cross product, mass 
properties, and EOM subroutines must be present. The model provides MR, 
A, and F as output. 


2.1.14.4 COORDINATE SYSTEMS 

The body coordinate system is used in this model. The RCS thrust axes are 
parallel to the body axes. Its X axis points toward the nose through the 
main propellant tank centerline, the Z axis points down in the engine gimbal 
pivot plane, and the Y axis points toward the pilot's right. 



2.1.15 


THRUST (Thrust Model) 


2.1.15.1 PROGRAM DESCRIPTION 

The inputs to the model are ambient atmospheric pressure, throttle 
setting, and engine gimbal angles and the outputs are forces and moments 
due to thrust from all engines in body coordinates. 


2.1.15.2 MATH MODEL 

ET i = PT i< T SLi + A 1 < PA - P » 


where: 



ET.j = engine thrust 

P-j... = throttle setting 

TsLi = sea level thrust 

A.j = engine area 

PA = sea level atmospheric pressure 

P = ambient atmospheric pressure 

i = engine index, maximum value = 12 


( 1 ) 


1 

*1 — 
X 
CD 
1— 

1 

r- 

TBY. 


l 

= 

Ltbz,J 



COS 0p 
0 


-si no 

P 


or: 


0 

1 

0 


-sinOj 

0 

cose 


P J 


cose. 


-sine. 


-sinGy 

COSOy 

0 



1 

1 — 
LlJ 

1 


0 


— 1 

o 

i 


( 2 ) 



T BX i 
T BY i 


ET^ cos 6p cos Gy 
-ET- sin Gy 
-ET^ sin 0p cos By 


(3) 
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2.1.15 THRUST (Continued) 


f tbx 

zT BXi 

- > 

f tby 

l7 BYi 

(4) 

f tbz 

zT BZi 



where 

T BXi’ T BYi 1 7 BZi = anci z components, respectively, of forces 

due to engine thrust for engine number i 
eY = yaw engine gimbal angle 
eP = pitch engine gimbal angle 

Ftbx s Fjby 1 F TBZ = ancl z components, respectively, of total thrust 

forces in body coordinates. 

MTXB. = TBZ. (ELY. - YCG) - TBY- (ELZ. ZCG) 

MTYB i = TBX. (ELZ i - ZCG) - TBZ. (ELX. - XCG) (5) 

MTZB i = TBY i (ELX i - XCG) - TBX i (ELY^. - YCG) 

MTXB^ , MTYB^ , MTZB^ = X, Y, and Z components, respectively, of moments 
due to engine number i. ELX^ , ELY , ELZ^ = X, Y, and Z components, 
respectively, of engine locations in body coordinates. XCG, YCG, ZCG = 

X, Y, and Z components of location of center of mass of the vehicle. 

MTB(l) = lMTXB. 

MTB(2) = zMTYB . (6) 

MTB(3) = EMTZB i 
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2.1.15 


THRUST (Continued) 


MTB(l), MTB(2), MTB ( 3 ) are X, Y, and Z components, respectively, of 
total moments due to engine thrust. 

All of the models described in equations (1) through (6) are in body 
coordinates. Body coordinates are fixed with respect to the vehicle. 

The X axis normally points forward along the main propellant tank cen- 
terline, the Z axis points down in a plane containing the engine gimbal 
pivots, and the Y axis points' toward the pilot's right forming a 
right-handed triad. 

2.1.15.3 INPUT/OUTPUT 

The parameters which must be supplied to the model as input are ELX^ , 

ELY^ , ELZ^ , PT i , T sl ^ , P, e p , 6y, XCG, YCG, and ZCG. The resulting 
quantities calculated by the model are FTB and MTB. The flight software 
command accepted by the model is PT. . The model has no calling arguments 
and calls no subroutine while operating. The Maspro, Aero, Autopilot, 
and Initialization subroutines provide the model inputs and the EOM 
accepts the model outputs. 
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2.1.16 TSHAPE (Trajectory Shaping) 

2.1.16.1 PROGRAM DESCRIPTION 

The trajectory and control parameters which must be calculated to accom- 
plish trajectory shaping are a D , desired angle-of-attack; e^, desired 
pitch attitude angle; 6^, desired engine deflection angle; and Z^g, 
desired body sensed acceleration. The values of e^, 6g, and Z DG g are 
dependent on a^. The boost flight is divided into three phases: vertical 
rise, tilt maneuver, and alpha policy. For each of these flight phases, 
exp is calculated differently. 

2.1.16.2 MATH MODEL 
2.1.16.2.1 GENERAL CALCULATIONS 


C Z a r 


= C^a 180 


q S C 


ZO 


N\ 


q S C 7 a 
Z t 


c (c m /c ) + x rr - x t 

mo zo CG F 


IP 


E (C m a/C z a) + X CG - X f 


'ig 


F ax Z CG 

X rr - X 
CG o 


'ig 


z rr - z n 

CG 0 i 

2 ? ** 

<% + z lg > 
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2.1.16.2.2.2 TILT MANEUVER (Parking Lot Tilt) 



The maneuver modeled here is a modification of the tilt maneuver pre- 
sented in the previous memo. Here, a m is used as the value of angle- 
of-attack at a time half-way between the time to begin and end the tilt 
maneuver, Jj and T^. Therefore, this procedure can fit any part of a 
parabola to the three points, depending on the value of a m , and its rel 
tion to ctg and c^. This is illustrated in Figure 2 t2. 
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2.1.16 



TSHAPE (Continued) 



where a m j, and are three examples of values for ^ which 


cause a different shape curve to be fitted. The calculation for 
during this maneuver is as follows: 

. 

calculated only at t = Tj 


l 0 


'OP 


tan 


'BZ 


'BX 


T + T 
= .T ‘D 


a D " a 0 


A. 


{a D - + (a D " a m ) 


T - T 


OP 


5 < t d - V 


« + A.Z + A,Z‘ 
ml 2 
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2.1.16 TSHAPE (Continued) 


2.1.16.2.2.3 ALPHA POLICY 


During this portion of flight a choice of several criterion for calculation 
of angle-of-attack can be made. 


2.1.16.2.2.3.1 


ZERO LIFT 
_C Z0 


C^a 


2.1.16.2.2.3.2 ZERO ANGLE-OF-ATTACK 


= 0. 


2.1.16.2.2.3.3 GRAVITY TURN 


For this criterion, gravity turn is defined as no’ acceleration normal to 
flight path. " 


'2B 


L B F B sin (6 


B 


6 CGB^ 


yy 


V 3B 


'IP 


— Sin <$„ 
m B 


-1 N 1 
VIP P 


yy 


'2P 


u F 

2 . 2v 2 h 0 


* (x lg' 'V 1 I 


yy 
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2.1.16 


TSHAPE (Continued) 


-N' p + F t 


IN - F (Z rr - Z„) 
oo ax v CG R' 


C 2P 6 CGO + 0 " C 2B 


Z o " K 3B 


B 1 K 3 ~ B 2 C 2P 
C 1P K 3 " C 2P K 2 


2.1.16.2.2.3.4 AERODYNAMIC MOMENT CONTROL 


^o N o + D z 
hp N ‘P 


2.1.16.2.3 REMAINING SHAPING AND CONTROL PARAMETERS 


The following parameters are calculated after a value for has been 
determined by some specified alpha policy. 
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2.1.16 TSHAPE (Continued) 


= tan 


1 1 ~ Vrz 

_ V RX_ 


= ci - e 
D R 


=1 o N o- F ax < Z CG ■ V 


= -1 N 1 
! 1P P 


(°o ' C 1P V ! yy + S CG0 " L B F B^B ‘ 6 CGB^ 


— N + N’ n a D - F sin 5 n - F D sin & D 
m o P o D B B 

_ — 


2.1.16.3 DEFINITION OF SYMBOLS 


Vari able 


B j j B 2 

C 1P’ C 2P 5 C 2B 


C , C 
mo m 

C f 
L ZO’ L Za 


Defi ni ti on 


Temp, variables 
Temp, variables 
Temp, variables 


Mean aerodynanic cord 


Aero, moment coefficients 

Aero, normal force coefficients (pitch) 


95 


2.1.16 TSHAPE (Continued) 


Variable 



K 2’ K 3’ K 3B 

L 0’ L B 
m 


N o 

N'p 

q 

s 

T 

t d 

T 0P 


V BX’ V BZ 
V RX’ V RZ 

Y v 7 

CG’ CG’ CG 
X, , Z, 

ig ig 

X o’ Z o 
X R’ Z R 


z 

Z DCG 


Defini tion 

Drag moment due to CG offset 
Total force acting on vehicle 
Aeoo. axial force 
Moments of inertia 
Temp, variables 

Moment arms from CG to orbiter and booster engines 
Mass 

Normal force at zero angle-of-attack 
Partial of normal 
Dynamic pressure 
Aerodynamic reference area 
Present time 

Time to end tilt maneuver (begin alpha policy) 

Mid point between T^ and lj 

Time to begin tilt maneuver (end vertical rise) 
Velocity components in body coordinates 
Inertial components of velocity relative to air 
Location of center of gravity 
Moment arm from CG to orbiter engines 
Average orbiter engine location 
Location of aerodynamic reference 
Location of body fixed accelerometer 
Temp, variable 

Desired body sensed acceleration 
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2.1.16 


TSHAPE (Continued) 


Vari able 


a. 


m 


6 CGB’ 6 CG0 


Definition 

Desired angle-of-attack 
Value of angle-of-attack at T^p 
Angle-of-attack at Ty 

Deflection for thrust through center of gravity 
at booster and orbiter engines 

Desired orbiter engine deflection angle 

Commanded pitch angle 

Acceleration at zero angle-of-attack 


2.1.16.4 INPUT-OUTPUT 

The inputs necessary for TSHAPE are the aerodynamic constants, configuration 
constants, I, m, CG's, T, Tp, Jj, Vg, V R , and a m - The outputs of this 
program are a Q , e c> 6 Q , Z DG , S GGB , and <5 CG q. 
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2.1.17 TVC (Thrust Vector Control) 

2.1.17.1 PROGRAM DESCRIPTION 

This model describes the motions of massless engines with limits on 
deflection, deflection rate, and acceleration. 


2.1. 

17.2 MATH MODEL 


a. 

5 NEW = 5 command 


b. 

’ 6 = ( 6 new ” ^old/ Atinie 


c. 

Test and limit 6 to obtain 

6 NEW 

d. 

6 = ^NEW " a old^ Atime 


e. 

Test and limit 6 to obtain 

6 new 

f. 

a new = a old + 6 new x Atime 


g- 

6 new = 6 0LD + a new x Atime 



6 NEW = Present en gioe deflection 


6qld = Previous engine deflection 


6 = Engine deflection rate 
6 = Engine deflection acceleration 

NOTE: 6q^p , 6ppp are only reset at the end of a computation cycle 

III. INPUT/OUTPUT 

Inputs: Pitch and yaw deflection commands from the flight control system 

(for each engine), deflection rate, and acceleration limits. 
Outputs: <S NEW , 6, 6 
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taa 


2 . 1.18 


INITIAL POSITION MATH MODEL 


2.1,18.1 PROGRAM DESCRIPTION 


This program calculates the difference between geodetic and geocentric 
latitude and uses it to calculate the initial state vector. This 
calculation needs to be done once each time either the launch 
azimuth or latitude is changed, 

2.1.18.2 MATH MODEL 


i|/L = Arctan \ (l-f ) 2 tan $ L 


R e ^" f) 

R, = r -J ~~\ + ho 

L v l-f ( 2 -f ) cos^ 


0 

R 

R 

R 

V 

V 

V 


X 

Y 

Z 

X 

Y 
Z 



= ‘f’L - <!^ l 

= R^ cos 0 
= R^ sin 0 sin A z 

= -R^ sin 0 cos A^ 

= 0 • ' - 

= a) R^ cos A^ cos 
= w R^ sin cos 
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2.1.18 


INITIAL POSITION MATH MODEL (Continued) 


Where: 


to = rotation rate of earth 
<j). = launch geodetic latitude 

f = earth flattening constant 
R g = earth equatorial radius 
A z = launch azimuth 

= launch geocentric latitude 
R l = magnitude of initial position vector 
3 = difference between geodetic & geocentric latitude 



initial position vector in platform coordinate 


initial velocity vector in platform coordinate 
altitude of vehicle CG above Fischer ellipse 


2.1.18.3 INPUT/OUTPUT 

The constants needed by this model are w, f, R g , and A^, The 
output of the program is R^,, Ry, R^, V^, V y, and V^. 

2.1.18.4 COORDINATE SYSTEM 

Platform System - an orthogonal system with its origin at the center of 
the earth, X axis parallel to the launch site gravity vector and positive 
in the direction opposite to gravitational acceleration. The Z axis 
lies in the launch plane and points downrange and the Y axis completes 
the right handed triad. 
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2.2 Rigid Body Boost Simulation (RIBBS) 

Between the time of the initial model development and when the SSFS first 
stage boost simulation became operational, boost dynamics and control analyses 
were performed using RIBBS. In the course of these analyses a number of 
improvements and modifications were incorporated into the RIBBS simulation. 
These improvements were subsequently incorporated in SSFS. Some of the 
improvements were trajectory shaping,, control gain calculation and state 
vector calculation subroutines. 

RIBBS was used extensively for analysis of the LOX-propane configuration 
which used liquid injection thrust vector control. For this study the thrust 
vector control scheme in RIBBS was revised to take out engine gimbals and 
add thrust vectoring by use of propane injection. The multi-stage guidance 
system described in Reference 15 was also initially originated and checked 
out in the RIBBS program. RIBBS was also modified to provide plotting and 
printout capabilities in case of any type of termination during a trajectory 
run, by means of a subroutine called RESTART. 


101 . 


2.3 


POINT TIME STABILITY ANALYSIS USING MDELTA 


The Boeing developed MDELTA program was used to perform point time 
stability analysis of pitch and lateral directional vehicle axes. The 
MDELTA program was modified to output Bode (magnitude and phase vs. 
frequency) plots and Nichols (magnitude vs. phase) plots. These plots 
are extremely helpful for quick analysis in determining stability 
characteristics of a system. 

2.3.1 PITCH PLANE DYNAMICS 

A derivation and discussion of the pitch axis dynamic model is found i 
Reference 16. Fig. 2- 3 shows the control diagram for the pitch axis, 
and the equations-of-motion are shown below. 

8 = c m a c + ( X CG " V C Z a I 9 ? ' f~ ^ X CG “ V 6 E 

y y 

q, c 7 a • . p 

Z = Z + ve - - 6c + (g sin ej e 

mV m E 3 o 

A test case analysis of pitch plane dynamics was run using data from 

the Shuttle Delta-Delta configuration. These values are shown in 
Table 2-1. For the three time points selected (0, 80, and 200 sec.). 
The resulting Bode and Nichols plots are shown in Figures 2-4 through 
2- 9 
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FIGURE 2-3' RIGID BODY PITCH PLANE MODEL 











TABLE 2-1 


MDELTA 


VARIABLE 

NAME 

0 SEC. 

80 SEC. 

200 SEC. 

X CG 

XCG 

144.73 

129.94 

103.54 

X R 

XR 

235. 

235. 

235. 

C Za 

CZA 

-.062 

-.0548 

-.0221 

c ma 

CMA 

-.1149 

-.1125 

-.04283 

c 

CBAR 

74. 

74. 

74. 

q 

QUE 

0. 

579.18 

9.71 

s 

AREA 

8308. 

8303. 

8308. 

'y 

IY 

482136000. 

439268000. 

309800000. 

u 

V 

V 

0. 

1470.735 

10099.187 

F 

F 

6480000. 

6996500. 

2315700. 

X E 

XE 

0, 

0. 

0. 

m 

MASS 

149948.1 

109467.5 

53119,9 

G 

G 

32,2118 

32.0787 

31.5253 

6 o 

THETO 

-.0354 

-.70404 

-1.679966 

K 0 

KTHET 

.5 

.2 

.4973 


KTDOT 

» 

.7 

.7612 

.6972 

k a 

KA 

0. 

.007 

0. 

% 

SIGN 

1. 

1. 

1. 

to 

n 

WN 

30. 

30. 

30. 

w c 

WC 

20. 

20. 

20. 


PITCH PLANE DYNAMICS DATA 
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ANGLE CDEGD MAGNITUDE CDBD 



RIGID BODY PITCH PLANE AT T = 0. 
(KTHET=. 5. KTDOT=. 7. KA=0) • 





RIGID BODY PITCH PLANE AT T=0 



FIGURE 2-5 
1 06 


ANGLE CDEG) MAGNITUDE CDB) 





RIGID BODY PITCH PLANE AT T=80 
( K I HE! =.2» KTDOT ", 76. KA= . 007) 




FIGURE 2-6 
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RIGID BODY PITCH PLANE AT T=80 



FIGURE 2-7 
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ANGLE CDEG3 MAGNITUDE CDB) 







MAGNITUDE CDB) 



RIGID BODY PITCH PLANE AT T=2Q0 



FIGURE 2-9 
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2.3.2 LATERAL DIRECTIONAL DYNAMICS 


The Lateral directional equations-of-motion and lateral control system 
were developed in Reference 17 but never used because of time allowance 
in scheduling of the tasks to be done under this contract. The Lateral- 
Directional control system is shown in Figure 2-10, and the vehicle 
dynamic equations are as follows: 


T6y = m V *8 - fiq C £ - (sin ct) mV P + (cosa)mV r - (gm) © 


T X e 2 + Z e 2 « - <C, ) 6 - 
Sqb g 


*xx ^ip b r— r + - r 

P + P + ir Z\T r Sqb 


sqb 


2V 


+ (C, ) 6 


T s| 4 v=( C n > » + C „P fv P+ 


xz 


zz 


Sqb y 


sqb sqb 


+ C nr 2 V r + C n 6 a 

6 a 


0 = (cos e ) P + (sin e ) r 


• COS d> 

tb = r 

y cos e 

0 
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2.4 


LINEAR TANGENT GUIDANCE LAW (LTG) DEVELOPMENT 

The explicit guidance law presented in this section (LTG) is very similar 
to the Saturn Iterative Guidance Mode (IGM) of Reference* 53 . Anyone 
familiar with the IGM can follow through the development term by term 
and see that LTG could be considered as an improved and somewhat simpli- 
fied IGM. It can also be shown that E Guidance, the Apollo lunar ascent 
guidance, is a special case of the generalized derivation presented in 
Section 2.4.1. The guidance law presented in Section 9.4.2 of Reference 54 
is a modified form of Apollo E Guidance. 

LTG is characterized by its utilization of the linear tangent steering 
law resulting from application of classical optimization theory. The 
derivation of LTG is presented in this Section. The generalized form 
of the steering commands is shown in Section 2.4.1 without defining the 
specific form of the terminal position and velocity vectors and the gravity 
model. (The major difference between IGM and E Guidance is the difference 
in the form of the terminal position and velocity vectors and the gravity 
model.) A simple second order gravity model is derived in Section 2.4.1- 
Implementation of the gravity model is incorporated in the equation flow 
chart os Section 2.1.7. In order to desmonstrate the flexibility of LTG to 
applications other than boost, it is also shown in Section 2.4.1 how to use the 
scheme to provide a throttle command for rendezvous. 

A computational flow of multi-stage guidance equations, as implemented in 
RIBBS (Boeing Rigid Body Boost Simulation Program), is shown in Section 2.1.7. 
Multi-stage guidance was implemented for two reasons: (1) to accommodate 

a constant acceleration phase as well as a constant thrust phase in the 
orbiter and ( 2 ) to enable active guidance turn on during booster flight. 

Linear tangent IGM type guidance was implemented in RIBBS rather than E 
Guidance because E Guidance, as implemented in the MSC/GCD Space Shuttle 
Functional Simulator, was derived for a one stage vehicle and, therefore, 
cannot tolerate a discontinuity in acceleration. The main simplifying 
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2.4 LINEAR TANGENT GUIDANCE LAW (LTG) DEVELOPMENT (Continued) 

assumption of E Guidance is that the effective gravity (gravity plus 
centrifugal acceleration) can be represented as the product of thrust 
acceleration and a linear function of burn time. Effective gravity is 
continuous with time and decreases with the square of the velocity. 
Therefore, the E Guidance effective gravity assumption is not very ac- 
curate for a constant acceleration burn, since this means that effective 
gravity is approximated by a linear function of time (whereas it is 
actually a higher order function). However, E Guidance is quite accu- 
rate for a one-stage constant thrust burn, since the acceleration has a 
higher order shape. LTG, as shown in Section 2.1.7, eliminates these 
restricting assumptions. 

2.4.1 DERIVATION OF GUIDANCE LAW 

We will consider vacuum flight equations of motion of a one-stage vehicle 
in developing the basic LTG equations. In a general derivation it is not 
necessary to specify the form of (F/m), the thrust acceleration (e.g., 

F/m can be constant or thrust and mass flow rate can be constant). It 
is only necessary to assume that F/m is a time function that can be 
integrated in closed form. Also, it will be seen that implementing 
multi-stage guidance is simply a matter of expanding the closed form 
integrals of F/m. It is also sufficient to assume that the gravity 
vector can be approximated by a time function that can be integrated in 
closed form. 


2.4.1 


DERIVATION OF GUIDANCE LAW (Continued) 

The following definitions are used in the analysis of this section 

R Total inertial acceleration vector 

(F/m) Thrust to mass ratio 

g Gravitational acceleration vector 

* ie 

Ap Unit vector of Euler-Lagrange multipliers defining thrust 
direction 

V Inertial velocity vector 

Vy Desired terminal velocity vector 

VgQ Velocity-to-go (to be gained by thrust) 

R Inertial radius vector 

Ry Desired terminal radius vector 

R gq Distance-to-go (to be gained by thrust) 

Tqq Time-to-go (burn time plus coast time) 

o c Commanded inertial pitch angle 

\jj c Commanded inertial yaw angle 

The following integrals are defined assuming that T^q is given. 

It is sufficient here to assume that closed form integrals exist 
for these quantities. The expressions for the integrals are given 
in Section 2.1.7. 

* In the following analysis is used to denote a unit vector. 
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2 . 4.1 


DERIVATION OF GUIDANCE LAW (Continued) 


Ao 

L - J (F/m)dt 

o 





(F/m) s ds 


dt 



2 . 4 . 1 . 1 Steering Commands. 

The vacuum equations of motion can be written according to the above 
definitions as: 

i) R = £ A F - 9- (i .e. g h -g) 


It is desirable to determine the thrust direction Up) time history that 
satisfies the terminal velocity (Vj) and radius ( R p ) and simultaneously 
minimizes burn time.- 


According to classical optimization theory the direction cosines of the 
unit thrust vector are defined by linear functions of time where the 
six coefficients are constant Euler-Lagrange multipliers (i.e., for a 
flat earth assumption or expressed in a spherical earth local horizontal 
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2. 4. 1.1, Steering Commands (Continued) 

coordinate system rotating with the vehicle)*. Also, the value of one 
of the six coefficients is arbitrary due to homogeneity of the equations. 
Therefore the reference thrust vector can be expressed as a unit vector, 
since three of the coefficients define the thrust vector and the re- 
maining three define the rate Of change. The solution for the coeffi- 
cients is made simple by defining the thrust vector as an expansion 
about a reference unit vector, of the form: 

2) Ap = A + A (t“K), where 

K is a constant .time about which the linear expansion (equation 2) 

A 

is made. A is a unit vector and A is the time rate of change of V, 

In general, the downrange component of terminal position (R-j.) is 
unspecified. It will later be shown that this leads to the fact 
that Tis normal to A, or expressed as a vector dot product, 


A 

3) A '* A = 0. 

The unit thrust vector can be written as: 

' A + A~ (t-K) 



1 + A . A (t-K; 

since A is a unit vector and A . T = 0. Equation 4 is the key to the 
simplicity and accuracy of this guidance algorithm. 

r * 

If o is defined as the angle between a f and A, then by definition of 

A /V * 

the vector dot product Ap . A = cos o, since both vectors are unit 
vectors. Making use of equations 3 and 4: 

* In application LTG is piece wise linear with time, somewhat approximating 
the Euler-Lagrange equations of a spherical earth in an inertial coordinate 
sys tern . 
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2. 4.1.1 Steering Commands (Continued) 

1 

5> Xf ' * " Vl + X 2 (t-K) 2 ' = C0S6 

9 o 9 * # o p p 

where a £- A . A. Therefore, A (t-K) = tan 0 since cos 0 = 

1 1 
Vsec^e Vl + tanks' 


It will be seen later that K is approximately equal .5TgQ. From 
equation 5 it is seen that when t = K^.5TgQ, cos 0 = 1, or 0 = 0. 
(" ~" denotes “is approximately equal to"). Experience has shown 
that it is a very good assumption in the derivation of the guidance 
equations to assume that: 


6) cos 0 ~ 


A/i + A 2 (t-K) 2 


1, 



The above approximation is compensated for in the manner that it 

affects the calculated value of Tqq and the integrals. The net effect 
of the assumption yields more accuracy than a small angle approximation 

Substituting equations 4 and 6 into equation 1 and rearranging, the 
approximate equation of motion becomes: 


7) £[ 

in L 


A + A 



R + g. 
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2.4.1 . 1 Steering Commands (Continued) 


Specifying the form of Ap in equation 7 reduces the guidance problem 
to determining values of A, a”, and K that. satisfy Vy and R^, specified 
terminal conditions. Integration of equation 7 according to the 
above definitions yields: 


8 ) 


L A + A 


(J - LK) = V T 


- V + AV 

g 



and 


9) S A + A (Q - SK) = R t - 


VT G0 + AR g 



The downrange component of terminal position (let us call this the Z 
component, or R-^) is unspecified in equation 9; therefore in com- 
ponent form there are six equations in eight unknov/n^A|, A3, 

• « • ■ 

Ap, A3, Rj3> and K) . The solution for Aand K is found by setting 
J - LK = 0 in equation 8, from which: 


10) K = J/L and La = V GG or 


11) A = V gq /L and 

12) L = ABS (V G0 ) = V G0 

Time-to-go is determined from equation 12 as shown in Section ,2.1.7. 
As mentioned earlier K *t.5T G q. This can be seen by inspecting 
equation 10 and 
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2.4.1 .1. Steering Commands (Continued) 


K = J/L = 



.5(F/m)T 2 G0 

(Wt G o 



(i.e. F/m = (F/.m) average) 

If the acceleration (F/m) is constant then K is exactly .5Tgg. 

Equation 11 shows that the average thrust direction is in the velocity 
to-go direction. 

An additional scalar equation is obtained by performing the vector 
/% 

dot product of A and equation 9, from which: 


13) x. R gq = S + A . A (Q-SK) and 

A • 

14) = (S ~ A j ^qq2 “ ^2 ^G02^^3 ^ ^ ^ (Q-SKj/A^ . 


The quantity A. R^g is unspecified since it contains which is 
unspecified; therefore, there should be no component of x in the a 
direction, which means that A is normal to A as expressed in equation 3 


3) 



Equation 3 could be called a transversal i ty condition in the classical 
optimization sense. 


2. 4. 1.1 Steering Commands (Continued) 

At 

From equation 5 it was shown that if A.x =0 
tan 2 0 = A 2 (t-K) 2 or 

tan 0 = ± x(t-K) = ±-(-xK+At)Sa + bt. 

(In the inertial coordinate system: 

^1 + A^(t-K) A + Bt 

tan x = : r , a bi-linear tangent 

» 3 + : 3 (t-K) ~ C + Dt 

a ; /\ 

law, however, in the coordinate system defined by x and R^q - x(x . Rgg), 
tan e - a + bt, a linear tangent law.) 

Applying equation 3', equations 13 and 14 become: 

13) X . = S and . 


14 ) R G03 ( S " X l R G01 " A 2 R G02^ A 3 


Rearranging the Z component of equation 9 


15) R t3 - R 3 Rqq 3 + v 3 t G q " aR 9 3' 



Equation 15 is used to compute the range angle-to-go in Section 2.2. 
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2. 4.1.i Steering Commands (Continued) 

A and R gq are now defined in equation 9, therefore the solution for 
X becomes 


16) 


rr_ R G0 
A Q - SK 





This completes the solution for K, x, and a. The thrust direction 
is defined by: 



A + X(aT-K) 

V 1 + X^AT-IO 2 ^ 


or 



Up = A + x(aT-K), where 


aT = "T - T q (T is present time and T q is time when K, 

>■> 

A, and. x were computed). The steering commands are: 


18) 0 c = .tan" 1 (Upj/Upj) - n/2 and 

19) = tan" 1 [u F2 /(U 2 F1 + U^)* 

A summary of the equations in the order of calculations is now given. 
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2. 4. 1.1. Steering Commands (Continued) 


V G0 i V T - V + 1V g 


V G0 - ABS <W 


x = W V G0 


R G0i R Ti R i " V i T G0 + AR gi ^ 


R G03 " x i R 601 " X 2 R G02^ A 3 


K = J/L 


A = (Rg 0 - SA)/ (Q - SK) 


AT = T - T 


U p = A + a(aT-K) 


e c = tan -1 (U F1 /U F3 ) - s/2 


*c tan [ U F2^ U FI + U F3^ 2 ] 
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2.4.1 .2 Gravity Model 

The Saturn guidance law. Iterative Guidance Mode (IGM), of Reference 1 
has demonstrated that a sophisticated gravity model is unnecessary. 

•In the IGM, gravity is approximated as having a constant, average 
direction and magnitude, I.e. 

% = - 5 Cij + 9 ) . where 

9 AV = average gravity vector, 
g-j- = terminal gravity vector, and 
g = present gravity vector. 

Then the velocity and position terms due to grayity are expressed as: 
__ /GO 

20) dV g - = g dt = g AV T G0 , and 

o 

2 

dt - .5 g AV T go . 


However, a more accurate, yet simple, gravity model has been imple- 
mented in this scheme (as shown in Section 2.2) in order to be con- 
sistent with the accuracy of the steering commands. Referring to 
equations 7, 8, and 9 of Section 2. 4. 1.1, it is seen that the steering 
parameters (a and x) are on the left hand side and the gravity terms 
are on the right hand side. Section 2. 4. 1.1 showed the solution of 
the left hand side of the equation. The purpose of this gravity model 
is to make the accuracy of the right side const-stent with that of 
the left. 


21) 


AR g = 


_ f 


/GO r* 


J 


gds 
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2. 4. 1.2 Gravity Model (Continued) 

The gravity magnitude is defined as the average value 


22 ) 




= . 5g 



.where 


g = present gravity magnitude, 

R = present radius magnitude, and 
Rp = desired terminal radius magnitude. 


At insertion R = Rp, therefore, g^ = g. 


The average downrange and crossrange components of position are 
determined by assuming second order expansions of downrange and 
crossrange position. This will be illustrated by considering the 
downrange component of position (Z), 


23) Z(t) =• Z Q + Zt + .5 Zt' 


The coordinate system in Section 2.2 is such that Zp = 0. It is 
assumed that: 


1 = 1 


h - z 


average T 


, where 


GO 


Z-j. is terminal value of Z. Equation 23 becomes: 

(z T - Z) t 2 

24) Z(t) = Zt + .5 — y 


'GO 


2.4.1 .2 Gravity Model (Continued) 


Spherical gravity is assumed in the form: 


25) g z = g^y or substituting equation 24 in equation 25: 


• • 


26) g 7 = % 
R AV 


Zt + .5 ( Z T_2 
r G0 


- z)r ’ 


, where 


R^y is average radius magnitude, to be defined later. The gravity 
velocity and position terms are defined by equations 20 and 21. 
Integrating equation 26 twice from 0 to Tgg yields: 


AZ„ = 


g AV r 


g R 


AV 


Z T - Z' 


(1/2) z T Z g +(l/6)(4— ) y3 


t 60 y -go 


and 


AZ 


9 R 


9AV i(l/6)z Tqq + 


AV t 


(l/24)( 


Z T- Z i T 4 

T G0 i G °J 


, or 


27) AZ 


J AV 


•• • 


g 


r {2Z + V -f 
k av 


GO and 


28) AZ 


g AV 


g R 


AV 


(3Z + Z T ) 


'GO 
T y 24 


Substituting Tgg into equation 24 yields: 


29) Z(Tgg) = Zy = .5 (Z + Zy) Tgg , the terminal value of Z. 

Multiplying equations 27 and 28 by the left side of 29, and dividing by 
the right side gives: 
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2 . 4 . 1.2 


Gravity Model (Continued) 


" 9 /\y 

AZ_ = AV 


■g R 


A V 


(2Z + Z T ) [ >5 (2 + 2 t ) T j 


and 


‘V 


M (3Z + 2 ) 
r + 24 

k AV 


.5 (2 + z T ) r G0 


By manipulation the above equations become: 


50) 42 - (g AV I G0 ) [ 2 t (F g + 1)23 ] /R^ and 


31) 42 n - (.5 g AV T^[ Z T (2Fg + D/e] /R flv> 


AV : 


where 


• • » 


32) F g = Z/(Z + Z T ) 


The crossrange or y components of these quantities can be obtained 
by inspection of equations 30, 31, and 32, since the derivation is 
analogous. In the y direction Fq (equation 32) is: 


f g = y/(y + y T ) 

Consider a local coordinate system where the local radius vector defines 
the X direction, Z is parallel to the desired orbit plane, and y is 
orthgonal . In this system y-j- = 0 and y^ = 0, therefore: 

F G = 1 

y T (F G + l)/3 = ( 2/3 )y y 
y T (2 F g + 1 )/ 6 = (l/2)y T , 

127 



2. 4. 1.2 Gravity Model (Continued) 


So the average y position for gravity velocity losses is 2/3 of the way from 
local position to the desired orbital plane and the average for 
distance losses is halfway. However, the X and Z axis of the guidance 
coordinate system implemented in Section 2.2 lie in the desired orbital 
plane. Therefore, the average y for gravity velocity losses is 1/3 
of the way from the desired orbital plane to local position. In the 
orbital coordinate system y^ =0, so y^ is used instead of yj, 
is used for the average radial distance. 

In summary, the gravity model is constructed by defining the average 
radius vector for velocity losses and the average radius vector for 
distance losses in terms of the terminal radius vector. As implemented 
in Section 2.1 .7, "tiie actual terminal radius vector is used instead of 
using approximations such as equation 29. Equation 29 is used only to 
derive the factors 1/3, (Fg, + l)/3, 1/2, and (2Fg + l)/6 as shown in 
Section 2.1.7. 

Inspection of Fg = Z/(Z + Z^) shows that Fg ranges from 0 to 1/2 as 
Z ranges from 0 to Z^. Therefore, the above factors range-as follows: 

1/3, crossrange velocity factory's constant. 

(Fg + l)/3, downrange velocity factor, ranges from 1/3 to 1/2, 

1/2, crossrange position factor, is constant. 

(2Fg + l)/6, downrange position factor, ranges from 1/6 to 1/3. 

So, at orbital insertion the factors are 1/3, 1/2, 1/2, and 1/3, 

If these factors were all constant and equal to 1/2, this gravity model 
would be almost identical to the gravity model of the Saturn IGM of 
Reference 1. 
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2. 4. 1.3 Throttle Command for Rendezvous 

Although there is presently no requirement for rendezvous capability 
in the SSV boost- to-orbit, this section is presented to illustrate the 
flexibility of this guidance scheme. 

In general, rendezvous implies that all components of the desired 
terminal position and velocity are specified. However, if rendezvous 
with a target vehicle is being performed, the position and velocity 
of the target are functions of the burn time of the vehicle performing 
rendezvous, and an iteration is involved. But for simplicity the example 
given here can be thought of as lunar descent to a hover point over a given 
landing site. In this case the desired terminal position is a given 
altitude over the landing site and the desired terminal velocity magni- 
tude is zero with respect to the landing site. 

We will assume that the vehicle propulsion system has a constant specific 
impulse (or exhaust velocity) and that when closed-loop throttling begins 
a constant thrust level is desired for achieving the terminal conditions. 
(It is even easier to compute a constant acceleration command). 

Equations 12 and 13 are used to compute the time-to-go and throttle 
command, i.e., 

12) L = V Q0 and 

13) S = A.R go . 

The quantities in the above equations are functions of T GG ; however, it 
will be assumed that the equations are recycled through until T GG stops 
changing (as in the range free mode shown in Section 2.1.7). The 
expressions for the integrals, L and S, are given in Section 2.1.7. 

One stage constant thrust is assumed. 
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2. 4. 1.3 Throttle Command for Rendezvous (Continued) 
33) L = V ex ln [ x/(t - T gq )J = V G0 


S L T G0 ' J " L T G0 ■ Lt + V ex t G 0 or 


34) S = L (T g0 - t ) +V ex T gq = X . R gq 


Define aR h a. R 


GO •* 


Rearranging equation 33: 


- v ro /V Px 

35) T gq = t(l-e G0 eX ) * 


Defining e' = 1-e 


-W V ex 


36) r G0 - xe' 


Substituting into equation 34: 

VgQ ( re' - t) + V gx Te' ■= aR, from which 

37) T = AR/[v G0 (e'-l) + V ex e'] . ; • 

Remembering that t = V /(F/m) it follows that the throttle command i 

38) (F/m) = V /t , where t is computed from equation 37. 

Time-to-go is now computed by substituting the above value of r into 
equation 35. The steering angles are computed as in Section 2. 4. 1.1. 
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2.4.2 COMPARISON OF IGM WITH LTG 

In order to compare IGM with LTG the vector derivation of LTG is con- 
verted to an angular derivation and the linear tangent law is introduced. 

2. 4. 2.1 Steering Parameters i ^ 


Referring to Reference 1, the form of the steering angles is 

x p = x p " K 1 + K 2 t = x p + K 2^ t_i< ) ( where K = 0/L) 

X y = X y - K 3 + K 4 t e J y + K 4 (t-K) 

or 

1) x p - x p = K 2 (t-K) and 

2) Xy ~ Xy = K 4 (t-k) 

however, the linear tangent form is: 

3) tan(xp - x p ) = K 2 (t-K) 

A) tan( Xy - x y ) = K 4 (t-K). 

Differentiation of Equations 3 and 4 yields 

5) x p sec 2 C Xp - Xp) = K 2 and 

6) * y ^d Z (x y - x y ) = K 4 . 

Inspection of Equations 3, 4, 5, and 6 shows that when t = k, Xp - 
x y = x y’ X P = ^2 ’ anc ^ x y ^4‘ 

The relationship between vector and angles is: 

■sin Xp cos Xy ' 

sinxy 

COSXpCOSXy. 


7) 


X F = 


x + T(t-K) 

[l + X 2 (t-K)*] 


21 % 
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X t 


2.4. 2.1 Steering Parameters (Continued) 


Differentiation of the above equation yields: 


8) a/u - a (t-K)u = 


Xp c osXp c osx y - x y sinx p sinx y 
*y cos *y 

L-XpSinx p cosXy - x y cosx p sin Xy . 

2/. „ \ 2- u 


where u s A + A(t-K) and u e ABS(u) = [l + Ar(t-Kr] 


Inspection of Equations 3, 4, 5, 6, 7 and 8 (for t = K) shows that: 
9) 

10) A = 


ID 

12) 

13) i = 

14) 


sinx cosy 
P 


sin x . 


cosx p cosXyJ 


and 


"K 2 cosx p cosx y - K 4 sinx p sinxy 
K 4 cos x y 

--K 2 si n x p cos x y - K 4 cosx p sinXy. 


From the above equations 


- - •? 2 2- 2 
A • A E A ^ COS^x + K 4 . 


The first order .approximation of LTG and IGM could be thought of as an 
average value approximation: 



[>/ 


1 + 


(K 


2 - 
cos x. 


K 4 )(t-K) 2 


dt = 1, 


which turns out better than a small angle approximation, by the way that 
it is compensated for in other calculations. 
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2 . 4 . 2.1 


Steering Parameters (Continued) 


The IGM steering parameters, Xp> Xy> Kj , K 2 , and K^, are obtained from 
Equations 9 through 13. It was shown in Section 2.4.1 .1 that 

15) x - Unit(Vg 0 ) and 


16) x = 


R G0 " Sx 


Q - SK 

From Equations 9, 10, 11 and 15 it follows that 


17) 

x" p = 

tan (Vq 0 i/ v g03^ 

and 

18) 

X 1 
<< 

11 

tan [ V G02 /( V G01 

+ V i 


From the definition K 3 = K^K and Equations 10, 13, 15 and 16 
(R GQ2 - S sinx y )K 


19) K 3 = 


cosx y (Q - SK) 


From the definition = K 2 K and Equations 9, 12, 15 and 16 


20) K 


R G01 “ S sin Xn C0S X 


£ ^ + K.tanx tanx ( 


cos Xp cos Xy (Q - SK) P y 


and 


K 2 - Kj/K 

A direct comparison can now be made between the steering parameters of 

LTG and IGM. The above equations for x n > x v > K,. and K. are identical to 

p y o 4 

those on pages 38 and 39 of Reference 1. However, the equation for K^ 

on page 40 omits the term K^tan^tanx of Equation 20. This equation 

o p y 

should be 

21) Kj = K x + K 3 tanx tanx . 
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2. 4. 2.1 Steering Parameters (Continued) 

The equations 

x p = x p " K ! + K 2 At and 

Xy = Xy -, K 3 + Vi 

should be 

22 ) x p = x p + tan _ 1 (-K 1 + K 2 At) and 

23) x y = x y + tan_1 (- K 3 + 


2. 4. 2. 2 I Range Angle 


An I GM range angle is implemented slightly differently than that in 
Section 2.2 due to the difference in the implementation of the guidance 
coordinate system. A local coordinate system is implemented in 
Section 2.2, whereas a terminal coordinate system is in Reference 1. 
Also opposite signs. are. on the gravity vector. An IGM, LTG type range 
angle calculation should be (employing the notation of Reference 53 ) : 

24) V 60 = (AX 2 + aY 2 + AZ 2 )^ 


25) 

26) 


7 = (c _ AX p _ AY p 

l G0 v V g0 k G01 V go G02 

AZ = Z G0 + ZT G0 + - 5g Rz T G0 
4> gg = tan _1 (AZ/X) 

4 . = tan _1 (Z ' 1 ' '/X 1 1 1 ' ) + <J> G0 


'GO 


AZ 


The range angle calculation on page 36 of Reference 53 is adequate for a 
first pass estimate. 
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2 . 4 . 2. 3 


Gravity Model 


A second order IGM gravity model is implemented differently than in 
Section 2.1.7 due to the difference in the coordinate systems. A second 
order gravity model is recommended for the sake of convergence and 
terminal accuracy (i.e., chi tilde steering is unnecessary). The nota- 
tion of Reference 1 is used in the following equations. Using the IGM 
coordinate system, the factors (F^ + l)/3 and (2F G + l)/6 of Section 2.1.7 
should be 1 - (F^ + l)/3 and 1 - (2F G + l)/6, where F fi = Z/(Z + Zy) . 

The equations are as follows: 

g = -(X 2 + Y 2 + 1 2 )\ 

y g g g pL 

Uav = ’^9 C 1 + J 


27) g DR = g AV 7 V x + (Y/3) 2 + [Z ( 2 


F g )/3] 


28) g. 


29) g vy 

30) g vz 


= 9 dr (V/3) 

= g n pZ(2 - F r )/3 


31) g DR = ** + (Y/2) 2 + [Z(5 - 2F fi )/6] 


32 ^ g RX g DR X 

33 ) g RY = g DR (v/2) 


34) g RZ = g DR Z(5 - 2F„)/6 


2. 4. 2. 4 Time-to-Go 


Time-to-go correction comes from the equation 


35) L = V G0 . 



2. 4. 2. 4 Time-to-Go (Continued) 


Define: <5V = Vg Q - L, 


A time-to-go correction is 

-6V/V 

6T n = t T^ " e 


exn» 


Assuming that 6V/V is small and using the first two terms of the 
3 exn 

exponential expansion, e~ x = 1 - x, the above equation becomes 


36) 6T n - x T 6V/» exn 


Tj = t - or' ^ exn / a max n depending on whether the insertion burn is 
constant thrust or constant acceleration. The present IGM time-to-go 
correction employs a square root algorithm 


. P . P . p 

w _ , /AX' + fiY + AZ ■, , x 
V™ ~ hi j- + L t ) 


'GO 


instead of 

• o . o •2 1 ' 
V G0 = + aY + AZ ' *• 
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2.4.3 COMPARISON OF E GUIDANCE WITH LTG 


The form of E Guidance that is discussed here is that presented in 
Section 9.4.2 of Reference 54. It is a modified fori\j of Apollo lunar 
ascent guidance. E Guidance can be derived as in Section 2. 4. 1.1. It 
is characterized by employing a local coordinate system rotating with 

O 

the vehicle, such that centrifugal acceleration (w R) appears in the 
equations of motion. The vector equation of motion is 

1) R = axr - - w 2 R)t, (spherical earth) 

r R 

where 

a = thrust acceleration 

A 

Xp = unit vector in thrust direction 
* 2 

g = -(-4)- - w R) = effective gravity, and 
R^ 

i = unit vector in radial direction. 

— I 

The simplifying approximations are 

A # 

2) Xp = u + ut (linear function of time) and 

3) g = a(g Q + gt) (product of thrust acceleration and a . 

linear function of time). 

The approximate vector equation of motion is 

4) R = a(u + ut) + (g + gt)^. 

o 

This equation can be written in the form 

5) r = a[x + x(t-K)] + ^-[g AV + g(t-K)] 

o 

where 

K = J/L = A 12 / A i i (Reference 2) 

% = (cj 0 + 9 K )lr’ and 
g = g^ 
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2.4.3 COMPARISON OF E GUIDANCE WITH LTG (Continued) 

The following definitions are made: 

s 0 = - U K + “o R o 

g p = -u/R 0 R D 

ap = final predicted value of thrust acceleration 

' _ 1 , a o v 

9 T~ (g F ' V 

l G0 h a F 0 

c = (q /a )i 1 

o ■ o' o'-V 

c = (g/a 0 )i r . and 
£ AV = E o + 

Equation 5 is now written in the form 


6) 

a [ a 

+ C AV 

+ {a + c ) ( t 

- K)] = 

R. 


Differentiation of 

Equation 6 as 

in Section 

2.1.1 yields 

7) 

L(X 

+ C AV^ 

+ (a + c)(J 

- LK) = 

*T - 

■ V = AV 

8) 

S(A 

+ C AV^ 

+ (A + c)(Q 

- SK) - 

r t - 

' R - v ' t go 5 AR 

Since K 

- J/L 

or J 

- LK = 0 




9) 

A * 

(a V - 

LCa V )/L 




10) 

^G0 

= aV - 

Lc AV’ 




11) 

V G0 

= ABS(V G q) , and 




12) 

L = 

V G0 
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2.4.3 COMPARISON OF E GUIDANCE WITH LTG (Continued) 

From Equations 7, 8 and 9 

13) A = (aR - SaV/L)/(Q - SK) - c s B - c 

Inspection of Equations 4 and 5 shows that 

14) u = a - xK and 

A _ Z. 

U = X 

from which 

15) u = ^ - c AV - BK + cK = - BK - c Q 

The third or downrange component of Ry is generally unspecified; however, 
in Section 2.1.1 it is shown that 

A • 

16) A-A = 0 and 

17) x-u = 1 - 

from which the third components of Equations 13 and 15 are derived: 

18) a 3 = -(AyAy + A 2 a 3 )/a 3 = u 3 and 

19) u 3 = (1 - AyUj - a 2 u 2 )/a 3 - 

The desired thrust vector is now defined as 

20) Uy ■ = u + u(t - t ) 

21 ) 3y - 3^Uy 

22) iyy = Unit (a T ) 

From Equations 17 and 21 it is seen that 

23) A • a T = a„ 

T o 


139 



2.4.3 COMPARISON OF E GUIDANCE WITH LTG (Continued) 


Equations 13 and 15 can now be converted to component form and compared 
with the equations in Section 9.4.2 of Reference* 54 . The following defi- 
nitions are made: 


a T 5 a o 


^eff = ^0 


A = aV^L - BjK, 

B = Bj 

c = av 2 /l - b 2 k. 

D = B 2 


These coefficients are identical to those in Reference 54. From the above 
coefficients and Equations 13, 15, 20 and 21, referring to the definitions 


eff 

= 9 0 and g: 


24) 

a Tr ‘ 3 t[ A + B(t '- V] 

- Vf - 9 (t " ‘o’ 

25) 

a Ty ' a T [ C + D <‘ - t 0 .] 


26) 

^GO ~ Vf - V - 

(Equation 46) 


A = Unit (Vqq) and 

27) a Jz = (a T - A^.^ - *2 a Ty^ A 3 (Equation 59). 
The desired unit thrust vector is 


28) 1 D = 


TR 


Ty 


Tz- 


1 


(a T 2 R + 3jy + £)* 


Equations 24 through 28 are in a linear . tangent form, whereas the equa- 
tions in Reference 54 are'inva linear sine form. The above equations are 
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2.4.3 COMPARISON OF E GUIDANCE WITH LTG (Continued) 

made identical to those in Reference 54 if the following approximations 
are made: 

29) g - 0 (Equation 24) 

30) A n c AV ='(i/2)T G0 £r e ff ■ (Equation 26) and 

31) a TZ = (a^ - a T 2 R - a 2 y )' 5 sign (Z 0 - Z) (Equation 27) 


Equation 30 is not a good approximation in the early part of boost flight 
A much better approximation is 


A 11 C AV = ^2a T ^ ABS ^t ‘ V). 9 eff • 


This approximation or a more exact form should be used in the time-to-go 
calculation.. Equation 27 yields more accurate, results than Equation 31 
and is simpler to implement since it eliminates logic. 

The form of E Guidance in Reference 54 and in this section is limited to 
one-stage guidance due to the gravity model. Effective gravity is ap- 
proximated as a product of thrust acceleration and a linear function of 
time; therefore, guidance cannot handle a discontinuity in thrust. Also, 
it is not a good approximation for a constant acceleration burn. A 
simple gravity model that eliminates these restrictions is shown in the 
next section. An alternate derivation of E Guidance is as follows: 

a = R + a [g AV + g(t - K) ] i^ 

\o - V T - V + Lg AV lR. K = J/L 

A = Unit (Vqq) , A • A ^ 0, S = A *Rqq 

^G0 = ^T ~ ^ V_T G0 + ^ Sg AV + g ^ ~ SK ^ h 
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2.4.3 COMPARISON OF E GUIDANCE WITH LTG (Continued) 

R G03 = ( S " A 1 R G01 " A 2 R G03^ X 3 
A = (R gq - Sa)/(Q - SK) , A 2 = A* A 

Ip = [a + A(t - t Q “ K)]/[l + ^ “ V 2 ] H 

Inspection of the above equations shows that the magnitude of A(t - t - K) 

A A . • U 

is the tangent of the angle between A and A + A(t - t - K) , and the 
magnitude of the latter quantity is the secant of the angle, since the 
magnitude of A is unity. 
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3.0 TASK II - RIGID BODY DYNAMICS AND CONTROL ANALYSIS 

Trajectory shaping techniques were developed as described in Reference 18. 
Drift minimum and loads minimum control laws were derived and continuously 
variable control gains were calculated using the methods reported in 
Reference 19. Design winds were synthesized using the techniques described 
in NASA Terrestrial Environment documents (1969 and 1971), References 20 
and 21, respectively. The differences between these references were found 
to be insignificant to the studies as discussed in Reference 22. The 
NR-GD 161B/B9T Delta Wing Booster/Delta Wing Orbiter and the MSC 036B/280 
inch SRM series burn configuration were analyzed during the first half of 
the contract period and reported in the semi-annual report (Reference 23). 
Since the semi-annual report, the MSC 040A/L0X propane liquid injectant 
TVC configuration and the parallel burn dual SRM configurations MSC 040 
and MSC 049 have been analyzed. These studies are summarized in the 
following paragraphs. 

161B/B9T 

A series of studies (reported in References 24 through 35) were performed 
to determine the control requirements of a Delta Wing Booster/Delta Wing 
Orbiter, typified by the NR-GD 161B/B9T. The unsymmetrical aerodynamic and 
mass characteristics of this type of piggyback configuration necessitated 
development of analytical techniques for trajectory shaping and control gain 
calculations. Because thrust vector control is inadequate to prevent large 
roll errors, several roll control techniques were investigated, including: 
limiting roll engine commands to permit large roll angles without loss of 
pitch and yaw control; aerodynamic roll control; and yawing to decrease the 
aerodynamic rolling moment. 

In the pitch plane, wind response studies revealed that drift minimum 
control can reduce trajectory dispersions about 65% as compared with a 
fixed gain attitude control law. The aerodynamic forces and moments 
were 10% higher using drift minimum but the force was still 1/3 less than 
the booster would experience during a 5g entry. The engine deflection 
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3.0 Continued 


requirement was approximately ±7.5 degrees, about half of which is due to 
winds and the other half required to control eg travel and no wind aero- 
dynamics. For crosswind control, ailerons were used during the high 
dynamic pressure region with thrust vector control only near lift-off and 
after 100 seconds of flight. Yawing away from the wind may also be useful 
for this configuration though the method of implementation would require 
additional study. Yawing into the wind tended to increase rather than 
decrease control requirements. 

036B/280 SRM 

The 036B/280 inch SRM booster series burn launch configuration was analyzed 
very briefly and found to be easily controlled in winds. Roll control used 
the orbi ter rudder. Reference 36 contains the results of this analysis. 

04QA/L0X Propane 

In February, 1972, the 040A/L0X propane liquid injectant TVC configuration 
was loaded into the SSFS Program as documented in Reference 37. A brief 
wind response study was performed, and it was determined (References 38 
and 39) that liquid injectant requirements could be reduced from 135,000 
pounds to 60,000 pounds by incorporating aerodynamic roll control. An 
additional study of load relief control was performed (Reference 40), and 
it was shown that additional injectant savings are possible by adjusting 
the slopes of the transition ramps from attitude control to load relief 
and back to attitude control. 
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3.1 SRM/040C CONFIGURATION ANALYSIS 


The SRM/040C configuration defined in Reference 41 was loaded into the SSFS, 
and both lift-off and inflight analyses were conducted. The lift-off 
investigation is described in Reference 42. The results of 135 flight 
simulations were reported together with the conclusions obtained from the 
investigation. Two modifications of the standard vehicle were considered 
using equations presented in Reference 43. The first involved tilting 
the SRM engines while the second translated and tilted the SRM engines. 

The results showed that all three configurations were flyable with reasonable 
values of thrust misalignment and unbalance. However, with an orbiter 
engine failure, only the configuration with the SRM's relocated to pass 
their thrust through the eg at lift-off was flyable during lift-off. 

The inflight analysis contained in Reference 44 included two parts: 1) con- 

trol law development including an active sideslip feedback roll control 
system; and 2) inflight studies where the vehicle's response to wind 
distrubances, SRM engine misalignments, thrust unbalance and engine failure 
were investigated. 

The control law development covered the following items: 

• Reviewing the available control moment authority in the maximum 
dynamic pressure regime. 

9 Selecting an engine deflection logic to maximize control authority 
without excessive roll -yaw coupling. 

e Modifying the basic SSFS control law to incorporate the use of 
rudder yaw control and rudder/aileron roll control . 

# Calculating attitude control gains 

9 Initial design of a roll command sideslip feedback control scheme. 
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3.1 Continued 

The inflight analysis produced several broad conclusions: 

0 "Simple" pitch plane load relief is inadequate for use with the 
sideslip control scheme - the control scheme itself minimizes qa 
and qB but positive attitude control is required to prevent 
excessive trajectory dispersions. 

0 The "back to the wind" sideslip control scheme was promising but 
it needed further development. 

0 The vehicle was satisfactorily controlled when flow in jimsphere 
winds. However, the more severe 95% synthetic design wind profiles 
sometimes produced unsatisfactory results. 

e The vehicle was not successfully flown with an engine out because 
dynamic pressure became excessive. An engine out guidance modifi- 
cation would be required. 
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3.2 SRM/049 PARALLEL BURN CONFIGURATION ANALYSIS 

In May, 1972, the SRM-049 configuration data was loaded into the SSFS 
program as reported in Reference 45. This vehicle was extensively analyzed 
both in lift-off phases of flight and in full boost to oribt dynamics and 
control studies as reported in References 45 through 51. In addition to 
the standard configuration with fixed SRM's, studies were initiated in 
July to investigate the use of vectorable SRM trim control. 

3.2.1 Lift-Off Analysis - Nonvectorable SRM's 

Results of investigations of the control capability and lift-off dynamic 
characteristics of the 049 launch configuration without vectorable SRM 
trim control have been documented in Reference 46. Lift-offs were flown 
in the presence of ground winds using various combinations of SRM angular 
misalignments and SRM thrust unbalances. Typical SRM thrust misalign- 
ments and the definition of the RMS roll yaw misalignment combination 
method is shown in Figure 3-1. 

Results indicated that the most demanding control problems occurred with 
SRM misalignments that produced vehicle perturbations of -roll and -yaw 
with a right crosswind. Plots of roll error, pitch error, and yaw error 
for the right crosswind case with 0.5 degree RMS SRM misalignments in 
-roll and -yaw are shown in Figure 3-2. 

Initially it was believed that with a right crosswind the most severe 
combination of roll yaw misalignments would be one with polarities of 
-roll/+yaw, which would be "in phase" with the wind disturbances. 

However, as demonstrated in Reference 46, the -roll/-yaw combination 
produces the more uncontrollable configuration. The reason for this is 
the "orientation" of the available TVC roll/yaw control torques. The 
-roll/-yaw combination places the control system nearer or outside the 
limits of control authority than does the -roll/+yaw combination. Likewise, 
due to symmetry, the +roll/-yaw combination is similar to the -roll/+yaw 
and the +roll/+yaw produces perturbations equivalent to the -roll/-yaw 
misalignment combination. 
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FIGURE 3-2 LIFTOFF RESULTS - 0.5° MISALIGNMENT - 
ATTITUDE ERRORS 
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3.2.1 Continued 


Results of simulations flown with 0.25° RMS SRM mi sa figments were, observed 
to produce no uncontrollable characteristics. Roll and yaw attitude 
errors for the .troll /-yaw and -roll/+yaw combinations of misalignments 
for all winds were found to be of the damped oscillatory form. Maximum 
engine deflections were on the order of 2 degrees, with maximum attitude 
errors of approximately 1 degree for these misalignments. 

For those simulations performed with -roll/-yaw combinations of 0.25° 

RMS SRM misalignments, it was observed that the engines would momentarily 
hit their limits in pitch, but as control was regained, engine deflec- 
tions were reduced. An example of these types of misalignments, attitude 
errors for the crosswind simulation with -roll/-yaw combinations of 0.25° 

RMS SRM misalignments are shown in Figure 3-3. As illustrated in this 
plot after approximately 6 seconds, control is regained and engine de- 
flections and attitude errors begin to decrease. 

Drift distances at tower clearnace for the simulations performed were 
found to be of acceptable magnitude, i.e., on the or^der of 2 to 3 feet. 
Exceptions to this were those flights with thrust unbalance; these cases 
produced crossrange drift distances approximately twice that of similar 
simulations without thrust unbalance. However, these were still felt to 
be within acceptable limits. 

3.2.2 Inflight Response - Nonvectorable SRM's 

Inflight response of the SRM/049 configuration to SRM misalignments, 
winds, and control system modifications has been documented in References 47 
and 48. The SRM-049 vehicle simulation was "flown" from lift-off to 
orbital insertion using a wind matrix composed of winds shearing to gusts 
at 10K, 20K, and 40K feet. Eastern test range 95 percentile winds were 
used as tailwinds and crosswinds; they reach a maximum velocity of 75 
meters per second. However, the maximum 95 percentile headwind velocity 
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3.2.2 Continued 


for the eastern test range is only 28 meters per second. Hence, the 42- 
meter per second Vandenberg Air Force Base wind was used for headwinds. 

In addition to the six wind velocity profiles investigated, consideration 
was given to inflight changes of wind direction or azimuth. Wind azimuth 
change profiles were prepared as directed by MSC/G&C Division. The wind 
matrix thus consisted of: fixed azimuth head, tail and crosswinds; winds 

that started as head, tail and crosswinds and then sheared away from these 
azimuths as the gust was approached; and winds that started an appropriate 
distance away from the nominal azimuth (both clockwise and counterclock- 
wise) and then sheared to the nominal direction at the gust. In all cases 
the wind azimuth remained constant both below and above the 2KM change 
region. Wind velocity and azimuth change profiles used are tabulated in 
Tables 3-1 through 3-7. 

Drift minimum control using accelerometer feedback as described in 
Reference 21 was employed in both pitch and yaw. The control gains were 
varied continuously to maintain a natural frequency of one radian per 
second and a damping ratio of one-half. This provided a resonant fre- 
quency of about 0.14 Hz which should be low enough to minimi-ze control 
interaction with structural bending and propellant slosh. No load relief 
methods were used in pitch or yaw in these studies. 

Four roll control techniques were investigated. The first was a simple 
attitude and attitude rate damping system. The second system added pro- 
portional and integral sideslip feedback to the roll attitude channel in 
order to roll the vehicle tail into the wind and decrease sideslip. The 
third technique was the same as the second except a. higher natural con- 
trol frequency was investigated. The fourth scheme was the same as the 
second except the proportional sideslip signal was fed to the roll rate 
channel. The fourth control scheme was selected as a new baseline for 
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HEADWIND PROFILES - I0K VANDENBERG, 
WESTERN TEST RANGE 
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HEADWIND PROFILES - 28K VANDENBERG 
WESTERN TEST RANGE 
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HEADWIND PROFILES - 40K VANDENBERG, 
WESTERN TEST RANGE 
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10K EASTERN TEST RANGE 
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CROSS AND TAILWIND PROFILES 
28K EASTERN TEST RANGE 
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3.2.2 Continued 

additional studies. These four control techniques and related parameters 
are summarized below. 





ROLL CONI 

fROL 

TECHNIQUES 




X 

Proportional 

Feedback 

Integral 

Feedback 

Proportional Feedback 
As Roll Rate Command 

1 

0.5 

0.5 

NO 


No 

No 

2 

0.5 

0.5 

Yes K =4 

Yes 

K = .4 T 
I 

No 

3 

1.0 

0.5 

Yes K =4 

Yes 

K .4 T 

I 

No 

4 

0.5 

0.5 

No 

Yes 

K = .4 T . 
■I 

Yes K =4 
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The baseline roll, yaw and pitch flight control techniques used are il- 
lustrated schematically in Figure 3-4. Notice that the rudder is used 
for roll control; it uses the same gains as for thrust vector control. 

The engines and the rudder were limited to +10 degrees deflection. 

With simple attitude control of roll (no sideslip feedback) all cases 
flew to orbit. The injected weight varied from 351,000 pounds for a 
headwind to 354,000 pounds for a tailwind (no wind is 352,500 pounds). 

Maximum dynamic pressure (q) and positive qa were observed for a constant 

O 

azimuth 40,000 ft headwind where q reached 569 Ibs/ft and qa reached 3848 
degree-lbs/ft . The time history of these variables is shown in 
Figures 3-5 and 3-6. (The machine plot miss the 3838 momentary peak.) 

The highest negative qa was -2964 (Figure 3-7) for the 28,000 foot constant 
azimuth tailwind. The engines just reach the 10° stops as shown in 
Figure 3-8. (The scalloped deflection curves reflect the linear inter- 
polation of the pitch command table.) Since the maximum headwind de- 
flection is only 8.5°, the engine bias could be shifted about 1° to -16°. 
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ORB ITER, PITCH; ENGINE DEFLECTIONS 





BETA - PSF 





TIME - SECONDS 
FIGURE 3-9 

INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - 
10K CROSSWIND - Q-BETA 
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3.2.2 Continued 


Maximum q 8 went as high as 4806 for the 28,000 ft crosswind that swings 

around to approach a tailwind. Roll control authority even with rudder 

2 

added is inadequate to hold qg above about 3000 degree-lb/ft . Roll 
control was lost during the high q region for all crosswind cases and 
the lower altitude azimuth shear tailwinds which approach crosswinds 
after the shear. Figures 3-9 and 3-10 are qg and roll error histories 
for the 10,000 foot shear crosswind. 

Proportional plus integral sideslip feedback (control scheme 2) was 
added, as shown in Figure 3-4, to roll the vehicle to a tail into the 
wind orientation. (The input to the integrator is opened to prevent 
control saturation if the engine deflection command exceeds eight degrees.) 
This modification decreased the 28,000 ft qg crosswind from 4806 to 910. 

The maximum qg for this scheme was 2980 for a 40,000 foot crosswind. 

Roll attitude errors remained as high as 54 degrees for control scheme 2 
because of slow vehicle response. The natural frequency was then in- 
creased to 1.0 radian per second (control scheme 3). The resultiwas to 
increase the roll attitude error (up to 83 degrees) because the scaled 
error signal reached the eight-degree limit and this combined with the 
increased rate gain limited the roll rate command capability to an even 
smaller value than previously. 

The proportional part of the sideslip signal was then transferred to the 
rate channel as shown in Figure 3-4 (scheme 4). This eliminated the 
sluggish response and reduced roll attitude errors to less than 10 
degrees . 

Table 3-8 summarizes the results of these studies. This table lists all 

2 

cases in which qa exceeded ±3,000 degree-lbs/ft . All cases were below 
4,000 and all were headwinds except for the 28K ft. crosswind without 
sideslip feedback. This case lost roll control in the high q region 
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3.2.2 Continued 


because of the high sideforce and turned belly to the wind similar to 
the 10,000 ft case shown in Figures 3-9 and 3-10. 

Table 3-8b lists all cases in which qp exceeds ±3,000 degree-1 bs/ft 2 . The 
principal cases were the crosswinds without sideslip feedback. With 
sideslip feedback none of the cases exceeded 3000 except for the 28K 
shear from a headwind case with the rate command scheme. For this wind 
profile there is no sideslip until the wind azimuth started to chage at 
about 22,000 feet. For this case the vehicle rolled through a high side- .. 
slip angle region. The similar lower altitude cases had a lower dynamic 
pressure so qp stayed lower, whereas at higher altitudes the wind velocity 
started to decrease soon enough to prevent $ from reaching a large value 
so qp remained below 3000. p also stayed lower for the cases without the 
rate command signal, although it is not clear why the slower response 
should have produced this result. 

Table 3-8c contains those cases where the weight at orbit insertion dif- 
fered from nominal, by more than 1,000 pounds. Of the twelve cases listed 
four had weight gains resulting from tailwinds adding energy. Seven of 
the eight weight loss cases where headwinds and the most. severe weight 
losses were headwind cases using sideslip feedback. Without sideslip 
feedback the headwind losses were less than 1500 pounds. . 

Table 3-8d lists the cases which had roll errors exceeding ten degrees. 
Large errors were experienced generally, until the proportional part of 
the sideslip signal was applied as a rate command. This control scheme 
improved vehicle response time sufficiently to keep the error less than 
ten degrees. 

SRM perturbations considered were thrust misalignments and thrust unbalance 
Lift-off studies had shown negligible problems with pitch misalignments 
of 0.25 degree so the inflight analysis was performed with combined 
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3.2.2 Continued 

roll -yaw misalignments of 0.1768 degree in each direction. (0.1768 x 2 = 

) 

0.25). The thrust unbalance consisted of 2.5 percent increased thrust for 
one SRM and 2.5 percent decreased thrust for the other one. 

The results are summarized in Table 3-9. None of the cases reached 4,000 
2 

degree-lbs/ft for q<* or qg and most of those exceeding 3,000 were head- 
winds with sideslip feedback. The only crosswind cases with q exceeding 
3,000 were those in which the misal ignments and thrust unbalance tend to 
force the vehicle to roll and yaw away from the wind. These cases are 
marginal and misalignments much greater than 0.25 degree would result in 
loss of control . 

Many cases showed weight losses exceeding 1,000 pounds. Only headwind 
cases exceeded 4,000 pounds weight loss and check cases of headwinds 
without sideslip feedback revealed no significant weight differences 
with and without roll -yaw misalignments. 

A gain of 4 was originally selected for the feedback signal. Parameter- 
ization of these gains is shown in Figure 3-11. From this graph it is 
seen that these values were well chosen. 

In summary, these studies demonstrated that, if SRM misalignments can 
be held to 0.25 degree or less and sideslip feedback is used when winds 
in the high q region are known to be more than 45 degrees away from 
headwinds, it appears that: 

2 

qa can be held below 4,000 degree-lbs/ft 

2 

qg can be held below 3,500 degree-lbs/ft 
weight to orbit can be held within 2500 lbs of nominal (for 
nominal thrust history). 
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TABLE 3-8 - CASES EXCEEDING CERTAIN VALUES - NOMINAL VEHICLE 
3-8a. qc< Greater than +3000 Degree-lb/ft^ 



No sideslip feedback oj - 0.5 

All headwinds - (Constant Azimuth are highest) 
1 OK Constant Azimuth Headwind 
28K Constant Azimuth Headwind 
40K Constant Azimuth Headwind 
28K Shear Clockwise from right crosswind 

Sideslip Feedback LO = 0.5 

Constant Azimuth Headwinds (same as above) 

40K Shear from headwind 
Sideslip Feedback = 1.0 

Constant Azimuth Headwinds (same as above) 

■ 40K Shear from Headwind 

Sideslip Feedback Plus Rate Command 1+7=0. 5 
Constant Azimuth Headwinds (Same as above) 

40K Shear from Headwind 
40K Shear to Headwind 

~8b. gg Greater than +3000 Degree-lb/ft^ 


No sideslip feed back- vJ =0.5 
1 OK Crosswind 
1 OK Shear from Tailwind 
28K Shear Clockwise from Right crosswind 
40K Shear Clockwise from Right Crosswind 
Sideslip Feedback = 0.5, 1.0 
None 

Sideslip Feedback Plus Rate Command 
28K Shear from Headwind 




3572 

3742 

3848 

3072 


3407 


3174 


3077 

3566 


M 


3598 

-3945 

4806 

4333 


3358 


171 


TABLE 3-8 (Continued) 

3-8c. Delta Weight to Orbit Greater than 1000 1b. 


No Sideslip Feedback 

1 OK Constant Azimuth Headwind -1480 lb 

1 OK Shear to Headwind -1480 

28K Constant Azimuth Headwind -1420 

28K Shear to Headwind -1420 

40K Constant Azimuth Tailwind +2300 

40K Shear from Tailwind +1230 

1 OK Constant Azimuth Tailwind +1030 

Sideslip Feedback = 0.5 


Constant Azimuth Winds (Same as above) 
Sideslip Feedback tO = 1.0 
Constant Azimuth Winds (Same as above) 


1 OK Shear from Tailwind „ -1480 

Sideslip Feedback + Rate Command i O = 0.5 
Constant Azimuth Winds (Same as above) 

1 OK Shear Clockwise from Right Crosswind +1285 

1 OK Shear to Headwind -4550 

28K Shear to Headwind -3540 

40K Shear to Headwind -1870 


R oll Error Greater than Ten Degrees 


No sideslip feedback CO =0.5 
All Crosswinds 62 -106° 

10K Shear from Tailwind 102° 

28K Shear from Tailwind 24° 

Sideslip Feedback OJ =0.5 

All Headwinds except Constant Azimuth 15 - 54° 

1 0 K Shear from Tailwind 19° 

Sideslip Feedback to = 1.0 

All Headwinds Except Constant Aximuth 14 - 83° 

all Crosswinds ' 16 - 43° 

1 OK Shear from Tailwind 67° 


Sideslip Feedback Plus Rate Command 6 0 = 0.5 
None 
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TABLE 3-9 - CASES EXCEEDING CERTAIN VALUES - SRM PERTURBATIONS (SIDESLIP 
FEEDBACK PLUS RATE COMMAND, 0.5) 

3-9a . qc< greater than + 3000 degree -lk>/ft 

28K Shear cw to headwind +Yaw, -Roll Misalign, +5% unbalance 3697 

40K Shear cw to Headwind +Yaw, -Roll Misalign, +5% unbalance 3808 

28K Shear ccw to Headwind - Yaw, -Roll misalign, -5 % unbalance 3566 

40K Shear ccw to Headwind. - Yaw, -Roll Misalign, -5% unbalance 3627 

(Headwinds without sideslip feedback /^50 higher than no misalignments) 

3-9b . tf-S greater than + 3000 degree -lb/ft^ 

28K Shear cw from headwind + Yaw, -roll misalign, +5% unbalance +3800 

40K Shear cw from headwind + Yaw, -Roll Misalign, +5 % Unbalance +3545 

28K Shear cw to Headwind -Yaw, -Roll Misalign, - 5 % Unbalance -3586 

40K Shear cw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -3243 

1 OK Shear ccw to Headwind - Yaw, -Roll Misalign, -5% Unbalance -3644 

28K Shear ccw to Crosswind -Yaw, -Roll Misalign, -5% Unbalance +31(34 

40K Shear ccw to Crosswind -Yaw, -Roll Misalign, -5% Unbalance +3194 

3-9c. Delta Height to Orbit Greater than 1000 lb 

10 K Shear cw from crosswind +Yaw, -Roll Misalign, +5% Unbalance +1320 
10K Shear cw, from Tailwind +Yaw, -Roll Misalign, +5% Unbalance -1125 
10K Shear cw to Headwind +Yaw, -Roll Misalign, +5% Unbalance -1160 

1 OK Shear cw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -5480 

10K Shear cw to Tailwind -Yaw, -Roll Misalign, -5% Unbalance +1095 

28K Shear cw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -4370 

28K Shear cw to Tailwind -Yaw, -Roll Misalign, -5% Unbalance +1060 

40K Shear cw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -3020 

40K Shear cw to Tailwind -Yaw, -Roll Misalign, -5% Unbalance +1.580 

1 0 K Shear cw to Left Crosswind -Yaw, -Roll Misalign, -5% Unbalance -1900 

28K Shear cw to Left Crosswind -Yaw, -Roll Misalign, -5% Unbalance -2030 

40K Shear cw to Left Crosswind -Yaw, -Roll Misalign, -5% Unbalance -1830 

10K Shear ccw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -3960 

10K Shear ccw to Tai Iwind -Yaw , -Roll Misalign, -5% Unbalance +1255 
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3-9c. (Continued) 

28K Shear ccw to Headwind -Yaw, -Roll Misalign, - 5 % Unbalance -2350 
28K Shear ccw to Tailwind -Yaw, -Roll Misalign, -5% Unbalance +1090 
40K Shear ccw to Headwind -Yaw, -Roll Misalign, -5% Unbalance -1250 
40K Shear ccw to Tailwind -Yaw, -Roll Misalign, ~ 5 % Unbalance +1510 
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3.2.3 Development of Vectorable SRM Control and Lift-off Analysis 

Vectorable SRM trim control was investigated in an effort to be able to 
trim out SRM angular misalignments on the order of 0.5 degree and 
larger. In an effort to rapidly implement vectorable SRM trim control 
into the 049 SSV dynamics and control simulation, the basic orbiter TVC 
baseline control system was utilized as a source of scaled attitude 
errors and attitude rates. Since the SRM's would be used only as trim 
control and would be of low frequency, it was assumed that the control 
gains associated with the orbiter TVC system would still be applicable. 

Lift-off boost dynamics and control studies on the 049 SSV utilizing 
several forms of SRM control logic, in addition to normal . orbi ter TVC, 
were performed as discussed in References 49 and 50. Five basic forms 
of control logic were investigated. 

Initial efforts to include SRM trim TVC capability into the Space Shuttle 
Functional Simulator (SSFS) employed use of integral attitude error feed- 
back. For each control axis the attitude error from the normal orbiter 
TVC system was scaled and integrated to provide trim control to the SRM's. 
Flight simulations were performed initially with no limits on the SRM 
integrators, but with one degree hardware limits on the SRM actuators. 
Flights were performed with 1 degree SRM misalignments in roll , with values 
of .02, .04 and .06 for the SRM integrator constant, and with values of .1, 
.3 and .5 for SRM actuator rate limits. Results indicated that loss of 
control occurred for the system due to "overcharging" of the SRM inte- 
grators. Since the actuators were limited to 1 degree and the SRM inte- 
grators were unlimited, this introduced excessive Tag into the control 
system. Limits were then applied to the integrators. However, large roll 
errors still persisted. 

Numerous intermediate control systems were investigated; however, final 
conclusions from the analysis (Reference 49) indicated that a control 
system using integral attitude error plus rate damping for roll and yaw and 
only integral attitude error in the pitch axis produced acceptable results. 
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3.2.3 Continued 


Figures 3-12 and 3-13 summarize the data obtained for those runs made with 
this system. Peak roll and yaw errors are seen to decrease as actuator 
rate limit is increased. In order to minimize power requirements to drive 
the SRM's, it was desired to keep actuator rate limits as low as possible 
and still hold peak roll to a reasonable value. From these plots a com- 
promise rate limit value of 0.3°/sec was chosen for a baseline system. 

This value along with an integrator constant of 0.02 were felt at this 
point to represent a satisfactory system. (An integrator constant of 0.04, 
although producing smaller roll errors, was found to lead to more "ringing" 
in the control system.) 

As directed by MSC/G&C Division, this system was modified slightly such 
that all control axes were composed of the integral of attitude error and 
rate damping for vectorable SRM trim control. This control system is 
illustrated in Figure 3-14. This addition to the control system presented 
a more unified approach to SRM trim control. Comparison of the results 
of this system with those of the previous one without rate damping in 
pitch indicate insignificant differences. 

Figures 3-15 through 3-18 present the results of simulated lift-offs of 
the 049 SSV configuration flown with 0.5° RMS -roll/-yaw combination SRM 
misalignments and ±2.5% SRM thrust unbalance. Shown in these plots are 
attitude errors, pitch and yaw orbiter engine traces, and SRM pitch gimbal 
traces. The same parameters are shown in Figures 3-19 through 3-22 for 
1.0° RMS -roll/-yaw combination SRM misalignments and ±2.5% thrust un- 
balance. 

3.2.4 Inflight Performance for Vectorable SRM System 

Inflight results were presented in Reference 51 verifying the capability 
of the vectorable SRM trim control system under severe wind conditions. 
Simulations were performed in the presence of 28,000 ft. crosswinds and 
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attitude errors (degrees) 


GIMB.SRM, .5DEG RMS C -ROLL-YAW INT=. 02*SRM RLIM= .3 



FIGURE 3-15 

LIFTOFF RESULTS - VECTORABLE SRM - 0.5° SRM 
MISALIGNMENT. - ATTITUDE ERRORS 
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Orbiter Engine Pitch Gimbal Traces 


GIMB.SRM* .5DEG RMSC-RGLL-YAW)> INT=. 02.SRM RLIM= .3 



FIGURE 3-16 

LIFTOFF PEStJLTS - VE CTORABLE SRM - -5° SRM 
MISALIGNMENT - ORBITER PITCH ENGINE DEFLECTIONS 
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Orbiter Engine Yaw Gimbal Traces 


GIMB. SRM» . 5DE6 RMSC -RGLL-/AW). INT=. 02»SRM RLIM= .3 



FIGURE 3-17 

LIFTOFF RESULTS - VECTORABLE SRM - -0.5° SRM 
MISALIGNMENT - ORBITER YAW ENGINE DEFLECTIONS 
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SRM Pitch Gimbal Traces 


GIMB. SRM» . 5DE6 RMSC -ROLL-YAW), INT=. 02* SRM RLIM= .3 



FIGURE 3-18 

LIFTOFF RESULTS - VECTORABLE SRM - .-0.5° SRM 
MISALIGNMENT - SRM PITCH ENGINE DEFLECTIONS 
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GIMB. SRM> 1 DEG RMSC -ROLL-YAW )> INT=. 02»SRM RLIM= .3 



FIGURE 3-19 

LIFTOFF RESULTS - VECTORABLE SRM -.1.0° SRM 
MISALIGNMENT - ATTITUDE ERRORS 
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Qrtri'ter Engine Pitch Gimbal Traces 


6rMB.SRM.i DE6 RMSC -ROLL-YAW), INT=. 02. SRM RLIM= .3 



FIGURE 3-20 

LIFTOFF RESULTS - VECT ORABLE SR M - 1.0° SRM 
MISALIGNMENT - ORBITER PITCH ENGINE DEFLECTIONS 
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Qrbi.ter Eng\ne Ginibal Trace? 


GIMB.SRM, 1 DEG RMSC -ROLL-YAW) , INT=, 02* SRM RLIM= .3 



FIGURE 3-21 

LIFTOFF RESULTS - VECTORABLE SRM 1.0 Q SRM 
MISALIGNMENT - ORBITER YAW ENGINE DEFLECTIONS 
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SRM Pitch Gimbal Traces 


GIMB. SRM, 1 DEG RMSC -ROLL-YAW), INT=. 02, SRM RLIM= .3 



FIGURE 3-22 

LIFTOFF RESULTS - VECTORABLE SRM 1.0° SRM 
MISALIGNMENT - SRM PITCH ENGINE DEFLECTIONS 
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FIGURE 3-23 - Roll Axis Channel of Control System 
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FIGURE 3-24 

INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - 
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FIGURE 3-25 

INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0 SRM MISALIGNMENT - Q-BETA 
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ORBITER ENGINE PITCH GIMBAL TRACES (DEG) 
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FIGURE 3-27 

INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - ORBITER PITCH 
ENGINE DEFLECTIONS 
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SRM PITCH GIMBAL TRACES (DEGREES) 





SRM YAW GIMBAL TRACES (DEGREES) 
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FIGURE 3-30 
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INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - RUDDER HORSEPOWER 
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: • FIGURE 3-32 


INFLIGHT RESULTS - NO SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - ACTUATOR HORSEPOWER 
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INFLIGHT RESULTS - SIDESLIP FEEDBACK VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - Q-ALPHA 
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FIGURE 3-34 



INFLIGHT RESULTS - SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - Q-BETA 
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FIGURE 3—35 

INFLIGHT RESULTS - SIDESLIP FEEDBACK - VECTORABLE 
SRM - 1.0° SRM MISALIGNMENT - ATTITUDE EERORS 
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INFLIGHT RESULTS - SIDESLIP FEEOBACK - VECTORABLE 
SRM 1.0° SRM MISALIGNMENT - ORBITER . PITCH ENGINE 
DEFLECTIONS 
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FIGURE 3-39 
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FIGURE 3-40 
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TABLE 3-10 


SIDESLIP FEEDBACK PARAMETERIZATION - 
28K CROSSWIND - VECTORABLE SRM 


28000 FT CROSSWIND 1° RMS ROLL YAW MISALIGMENT 

* 2° SRM DEFLECTION LIMITS - 2.5% THRUST UNBALANCE 


SIDESLIP 

FEEDBACK 

GAIN 

• 

q<* 

(PSF°) 

(PSF° ) 

ROLL 

- ATTITUDE 
(DEG.) 

WEIGHT 
TO ■ 
ORBIT 

(POUNDS) 

RUDDER 

HORSEPOWER 

HOURS 

ACT • TOR 
HORSEi TWER 
HOU S 

0 

2432 

4943 

-31.0 

352733 

.0064 


.3767 

.4 

-1381 

3582 

78.5 

352179 

.0034 


.3669 

.6 

-1980 

2683 

85.5 

352070 

.0027 


.3702 

.8 

-2287 

2130 

88.0 

351936 

.0025 


.3834 

1.0 

-2382 

1824 

88.7- 

351873 

• 

.0028 


.4042 


TABLE 3-11 

SIDESLIP FEEDBACK PARAMETERIZATION - 
28K HEADWIND - VECTORABLE SRM 

28000 FT HEADWIND 1° RMS ROLL YAW MISALIGMENT 

±2° SRM DEFLECTION LIMITS - 2.5 % THRUST UNBALANCE 


SIDESLIP 

FEEDBACK 

GAIN 

qa 

(PSF°) 

q/3 

(PSF°) 

ROLL 

ATTITUDE 

(DEG.) 

WEIGHT 

TO 

ORBIT 

(POUNDS) 

RUDDER 

HORSEPOWER 

HOURS 

ACTUATOR 

HORSEPOWER 

HOURS 

0 

3845 

-579 

-11.5 

345570 

.0004 

.3153 

.4 

3774 

-1208 

-34.1 

344950 

.0015 

.3404 

.6 

3705 

-1781 

-69.8 

344276 

.0028 

.3656 

.8 

3591 

-2557 

-119.0 

343648 

.0049 

.3959 

1.0 

3401 

-3229 

-161 .0 

343516 

.0070 

.4267 
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3.2.4 Continued 

headwinds because they had previously (Reference 48) been found to produce 
the most severe control requirements. 

The roll axis channel of the baseline control system was modified slightly 
for inflight analysis to allow for limiting of the sideslip feedback 
signal. This modification is shown in Figure 3-23. Flights were simu- 
lated without sideslip feedback and with sideslip feedback gains ( Kj 
of 0.4, 0.6, 0.8, and 1.0 with a beta signal limit (BUM) of 5°. 

Table 3-10 summarizes the crosswind results and Table 3-11 summarizes 
the headwind results. Several variables are shown, but the most interest- 
ing was qg. It is seen that, for the crosswind, qg decreases as the feed- 
back gain increases, while the opposite is true for the headwind. Based 
on these results, a feedback gain of 0.8 was selected as a baseline for 
future analysis. 

V 

Figures 3-24 through 3-41 show results of SRM gimbaling for the 28,000 ft 
crosswind simulations with 1 degree RMS misalignments. Figures 3-24 
through 3-32 are without sideslip feedback,, and Figures 3-33 through 3-41 
are with sideslip feedback. They illustrate the smooth response using 
SRM gimbaling and the small demand on the orbiter main engine gimbaling 
in the high q region. Comparison with Reference 48 showed that, for the 
headwind case with sideslip feedback, qg was reduced from 3566 PSF° to 
2557 PSF° when SRM gimbaling is added. For the crosswind case qg was 
reduced from 3194 PSF° to 2130 PSF°. 
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3.3 BASELINE CONTROL SYSTEM 


The baseline launch configuration control system software described in this 
section satisfies the requirement to control the SRM parallel burn vehicle 
from lift-off through SRM engine cutoff and orbital insertion. This con- 
trol system has evolved from the investigations of References 44, 48, 49, 
and 51. It is the result of efforts to provide active roll control for 
load alleviation and the option to vector SRM thrust such as to trim out 
SRM angular misal igments . Specific details presented in this section are 
for a nominal rigid body version of the MSC SRM-049 launch configuration. 
Compensation for bending and malfunctions can easily be incorporated. 

The system is adaptable to various control laws, including attitude con- 
trol (as in Saturn V), drift minimum control (with and without acceleraometer 
feedback), and/or load relief control. It also contains the option of 
using thrust vector control (TVC) on the SRM's as trim control to remove 
initial thrust., misalignments. In-addition, sideslip feedback may be em- 
ployed during boost phases of flight to roll the vehicle into a tail to the 
wind attitude as a means of load alleviation. 

The following paragraphs contain a description of the control systems 
with functional block diagrams of the pitch, yaw, and roll channels and 
a discussion of the configuration dependent variables. 

3.3.1 Software Description 

The baseline launch configuration control system is essentially a conven- 
tional booster control system. It is a digital autopilot receiving gimbal 
angles from a stable platform, body rates from body mounted rate gyros, and 
body Y and Z axes translational accelerations from body mounted accelerome- 
ters. It also receives prestored engine deflection and acceleration com- 
mands, and attitude commands either prestored or from a guidance system. 

From these inputs the control system determines and outputs deflection 
commands for the orbiter, SRM, and rudder actuators. 
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3.3.1 Continued 

The prestored commands for attitude, acceleration, and orbiter engine 
deflections are the measured values taken from a reference trajectory. 

The baseline reference trajectory was shaped for gravity turn (zero 
resultant acceleration normal to the longitudinal stability axis) in 
order that roll maneuvers might be accomplished without deviation from 
the desired trajectory. Variable control gains are calculated to satisfy 
the selected control Taw and desired rigid body frequency and damping 
ratio. For current studies the drift minimum control law was employed 
with a natural frequency of one radian per second with a 0.7 damping ratio 
in both pitch and yaw. Roll control utilizes both thrust vector and rudder 
control at a natural frequency of 0.5 radian per second and a damping 
ratio of 0.5. It should be noted that gain calculations were made assuming 
no SRM thrust vector control. However, since the SRM's are used for trim 
control (low frequency), this assumption is valid for the rigid body 
analysis. Development of gain calculation equations is given in 
Reference 19. 

A more detailed functional description of each channel of the control 
system is provided in the following paragraphs. Separate block diagrams 
for roll, pitch, and yaw channels are presented. 

3.3.2 Roll Channel Description 

Referring to Figure 3-42, the block in the lower right hand corner labeled 
"Vehicle and Subsystem Dynamics" represents actuator and sensor dynamics, 
prefiltering for the sensors, flight computer input and output registers 
and any A/D or D/A converters. Inputs to the control system from "Vehicle 
and Subsystem Dynamics" are the inertial platform outer gimbal angle and 
rate measurements from one or more body (or wing) mounted rate gyros. 
Outputs from the control system to the vehicle are orbiter engine, SRM 
engine, and rudder deflection commands. 
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ROLL CHANNEL - NEW BASELINE CONTROL SYSTEM 



















3.3.2 Continued 

Following the attitude loop in Figure 3-42, the platform gimbal angle (<j> ) 
is differenced with the inertial attitude command (<j>C) to produce the 
inertial attitude error (<j>E) which is transformed to a body roll error 
UEB) by adding the inner gimbal angle error times the sine of the middle 
gimbal angle (<|>EB = <}>EB + <f>EI si n^) . At this point the body roll error 
signal is differenced with the signal from the sideslip feedback inte- 
grator. 

The purpose of the sideslip feedback portion of the control system is to 
provide a means of alleviating large sideloads (qB) resulting from cross 
winds, by sensing the sideslip angle and producing appropriate control 
signals such as to roll the vehicle into a tail to the wind attitude. 
Referring to Figure 3-42, this is achieved by integrating the sideslip 
signal and subtracting it from the body roll error signal. In addition, 
the sideslip signal is scaled and^differenced with the roll rate signal. 
It should be noted that the sideslip signal is passed through a 0.2° 
deadband and limited. Following this is a block which represents a 
"ramp-on" and "ramp-off" scheduling. This is to allow the sideslip roll 
to be initiated after tower clearance and be turned off after the wind 
effect has become negligible (after the high dynamic pressure region). 

The switch prior to the sideslip integrator is activated by the limited 
orbiter engine roll commands. When the roll command becomes greater than 
8 degrees, the switch to the integrator is opened. This prevents the 
sideslip integrator from saturating the roll commands ana reserves 2 
degrees of roll control for rate damping commands. The body roll error 
signal, having been differenced with the sideslip integrator, is at this 
point split to supply commands to both the aerodynamic control and the 
thrust vector control. Referring to the roll channel block diagram, it 
should be noted that a limiter has been included in the attitude error 
loop. This limiter reserves some TVC authority for rate damping (even 
though the attitude error signal may be saturated) thus preventing 
oscillation and complete loss of control because of large signal in- 
stability. 


213 


3,3.2 Continued 

Returning to the "Vehicle and Subsystem Dynamics" block in Figure 3-42, 
roll rate gyro signals are filtered and blended as required. The roll rate 
error is then scaled and filtered separately for the engine and aero 
channels and is mixed with the attitude error signals to produce total 
error signals. The rudder deflection signal is limited to avoid over- 
driving the servos and sent to the output registers of the rudder. In 
the engine TVC loop the combined roll signal ( <5C) is at this point 
split to provide signals for orbiter engine control and SRM trim con- 
trol . 

For SRM trim control, the combined roll signal is rate limited, integrated, 
and sent to the SRM deflection logic. The limiter associated with the SRM 
integrator limits the integrator and its output. As noted in Reference 49, 
it was found to be necessary to limit the integrator rather than just its 
output to prevent it from "overcharging" and inducing excessive lag into 
the system. The baseline SRM deflection logic is simple, equal and 
opposite SRM pitch actuator commands. Finally, the SRM engine deflection 
signals are limited to the hardware limits and stored in output registers. 

In the orbiter engine TVC loop the combined roll signal ( <5C) is limited 
to the hardware limits and goes to the orbiter engine deflection logic. 

The orbiter engine deflection logic is complicated by the need to decouple 
yaw reaction from the roll command. Roll commands are sent to the appro- 
priate orbiter engines actuators to produce the desired vehicle torque, 
i.e., equal and opposite pitch actuator deflections to the lower two 
orbiter engines and a yaw actuator command to the upper engine. In order 
to remove the yaw induced by the upper orbiter engine, the lower engines 
are commanded in yaw in a direction opposite to the upper engine and by 
an amount equal to one half of the yaw deflection of the upper engine. 

This deflection logic is illustrated in. Figure 3-43. 
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3.3.3 Yaw Channel Description 

The yaw channel block diagram in Figure 3-44 is laid out in the same 
manner as the roll channel with the "Vehicle and Subsystem Dynamics" 
block in the lower right hand corner, the "Inertial Attitude Command" 
in the upper left hand corner and the "Engine Deflection Logic" on the 
right. The aerodynamic control circuit is missing because there is 
adequate engine control authority to handle yaw plane torques. 

Three features which were not included in the roll channel appear in the 
yaw channel. An accelerometer input has been added to the attitude 
loop, a forward loop integrator follows the limiter, and a rate command 
signal is differenced with the input from the rate gyros. 

The platform gimbal angle (ip) is differenced with the inertial attitude 
command (ipC) , transformed to body attitude error (ipEB = ipE I cos<t> - 
eEI sine}) cosi|>) and scaled. The error signal is then modified by summing 
an acceleration error signal. The acceleration error is obtained by 
blending appropriately filtered signals from one or more body mounted 
accelerometers and then applying gain and filtering compensation accord- 
ing to the control law or bending requirements. For the baseline rigid 
body design, the gain is varied with time to meet the drift minimum 
requirements of Reference 19. After provision for additional filtering 
(not used in the rigid body baseline) and after limiting the combined 
error signal, the forward loop integrator trims out steady state errors 
resulting from biases or off-nominal conditions. 

The rate command signal is introduced because the normal attitude command 
changes at a predictable rate. The rate command signal minimizes the 
tendency for the vehicle attitude to lag behind the attitude command due 
to rate feedback. If the vehicle rolls appreciably, the rate command 
must be applied partially to the yaw channel. In operation, the attitude 
command is differentiated, multiplied by the sine of <j>, scaled the same 
as the rate feedback and differenced with the rate feedback signal. The 
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YAW CHANNEL - NEW BASELINE CONTROL SYSTEM 













3.3.3 Continued 

rate error signal is then added to the attitude error signal and sent to 
the "orbiter engine deflection logic," where the lower two orbiter engines 
are commanded equally in yaw and also to the SRM trim control section. The 
SRM trim in yaw is the same as that previously described for roll. The 
command signal is rate limited, integrated and limited, and sent to the 
SRM deflection logic, where both SRM yaw actuators are commanded equally. 

3.3.4 Pitch Channel Description 

Comparison of Figure 3-45, which is the pitch channel block diagram, with 
Figure 3-44 (yaw channel) reveals that the only difference is the addition 
of prestored acceleration commands and prestored engine deflection commands 
Both of these command schedules are the nominal values from a reference 
trajectory. 

The prestored acceleration command is necessary because of the lack of 
vehicle symmetry, both in an offset center of gravity and in aerodynamic 
characteristics. The. results of these asymmetries is that the nominal 
normal acceleration (Z)'is typically non-zero. The addition of the pre- 
stored acceleration command permits use of accelerometer feedback with any 
feasible trajectory to modify the attitude error signal to conform to the 
desired control law (drift minimum for the baseline case). 

The prestored engine deflection command schedule provides coarse control 
of the engine gimbal trim in the pitch plane. Center of gravity travel 
alone could be handled by the integrator. Aerodynamic pitching moment 
increases sharply in the vicinity of Mach one and decreases again there- 
after. A high integrator gain would be required to follow these rapidly 
changing deflection requirements. Since a high integrator gain is un- 
desirable in the high dynamic pressure region, the prestored commands 
are used to control the nominal part of this variation. For the base- 
line case the integrator gain (t) is zero. 
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FIGURE 3-45 

PITCH CHANNEL - NEW BASELINE CONTROL SYSTEM 














3.3.4 Continued 

For the orbiter engine deflection logic the pitch channel is the same as 
the yaw channel with the exception that all orbiter engines are deflected 
equally in pitch. For the SRM trim deflection logic the basic principles 
are the same as presented before for the roll and yaw channels. However, 
the pitch SRM deflection logic is slightly different. Since ample pitch 
control is available from the orbiter TVC system but limited orbiter TVC 
roll control is present, the SRM deflection logic has been devised such 
that roll control takes precedence over pitch commands. If the combined 
roll and pitch commands to either SRM produce a command signal rate 
greater than the limit, then the rate to the pitch integrator is the rate' 
limit minus the roll rate command. Also, if the combined roll and pitch 
engine deflection command from the SRM integrators is greater than the 
limit, then the pitch deflection command is set equal to the limit minus 
the roll command. Thus full roll commands are preserved, and, if a 
reduction is necessary, it is taken from the pitch channel. 

3.3.5 Configuration Dependency 

The basic control system is designed for a dual SRM parallel burn con- 
figuration. The control gain and command schedules are, of course, 
configuration dependent. They will be affected by modifications to 
either the vehicle configuration or the reference trajectory. This will 
not impact the design, however, because all control gains and commands 
are continuously variable input functions of time. Some of them may 
subsequently be modified to functions of other trajectory variables, 
such as altitude and velocity, but this will involve only minor software 
modifications. 

Filters have not been developed. They will be functions of the vibra- 
tion and slosh characteristics of the launch configuration. The provision 
for incorporating filters at various points is intended to reduce the 
order of each filter. 
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4.0 TASK III FLEXIBLE BODY STABILITY AND CONTROL ANALYSIS 

Frequent changes in vehicle configuration prevented performance of 
flexible body analysis. However mathematical models were developed 
for use in the SSFS. These models were programmed for a parallel burn 
solid rocket motor configuration such as the 049 configuration analyzed 
in Section 3. The flexible body version of SSFS is currently being 
checked out. The models as contained in Reference 52 are presented 
in 'the following paragraphs. 

4.1 FLEXIBLE BODY PROGRAM DESCRIPTION 

This program contains the bending and slosh models for the launch 
configuration during first stage boost. It uses a generalized modal 
approach to bending which represents the elastic response by stan- 
dard normal modal equations with viscous damping. Included are models, 
for aerodynamic forces and moment^ and thrust forces and moments to 
account for bending effects as well as the tail wags dog contribu- 
tion to bending. The rigid body and elastic response equations 
provided here are uncoupled and are considered separately since the 
magnitude of the coupling is insignificant. The number of equations 
is very sensitive to the vehicle configuration and to the complete- 
ness of the bending analysis. Therefore, when data becomes available 
it is likely that only a small percentage of the general set provided 
here will actually be required for SSV analysis. 

The model sums all the forces acting on each of the equivalent mass 
points and for a given mode numerically integrates the sum with a 
second order linear differential equation in modal displacement. 

The number of mass points at which aero forces and modal displacements 

are calculated will be less than 50. The number of modes at these 

points will be less than 10 each. The number of slosh masses will 

be less than 5 and the number of modes per slosh mass will be less 

than 5. The EOM, guidance atmosphere and control subroutine must be 

present to provide inputs for this model. 
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4.2 Vibration Equations 


N1 


E (F . 6 . + F . d) . . + F -d)..) 

j _ 2 v axj 9 xij ayj azj zij 


+ Ml 

E (F. . <t> . . + F . • < 1 > . . + F, . d> . 
j = 1 tx J X1 J tyj 9 yij tzj 9 zij j 


+ K1 

E (F .6.. + F . d> . . + F .4..) 

■ _ j sx j xi j syj yij szj 9 zij ; 


+ N2 

E (M M . d.' . . + M .<*'.•) 

j = j \, axj 9 xij ayj 9 yij azj 9 zij ; 


+ M2 

E (M. . <b'.. + M. . d>'.. + M. . (J)'..) 
j _ 2 tx J x1 j tyj yij tzj 9 zij y 


+ K2 i i i 

E (M . (j) . . + M . d>„. . + M . d> . .) 

. _ 2 sx j xi j syj yu sz j z u 


= m. (q. + 2 a . q. + q.) 


q. = fq. dt + q. 
q. = / q. dt + q. 
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4.2 Vibration Equations (Continued) 


Where: 


N1 

N2 

Ml 

M2 

K1 

K2 

S* 


uj . 


m. 


axj 


ayj 

azj 

txj 

tyj 

tzj 


sxj 


syj 


szj 


M 


axj 


M 


ayj 


M 


azj 


= number of aero stations for aero forces 
= number of aero stations for aero moments 
= number of engines producing' forces 
= number of engines producing moments 
= number of slosh stations for slosh forces 
= number of slosh stations for slosh moments 
= modal displacement due to bending mode i 
= damping coefficient for mode i 
= frequency of mode i 
= normalized mass for mode i 
= aero forces in X direction at station j 
= aero forces in Y direction at station j 
= aero forces in Z direction at station J 
= thrust forces in X direction for engine j 
= thrust forces in Y direction for engine j 
= thrust forces in Z direction for engine j 
= slosh forces in X direction at station j 
= slosh forces in Y direction at station j 
= slosh forces in Z direction at station j 
= aero moments about X axis at station j 
= aero moments about y axis at station j 
= aero moments about z axis at station j 
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4.2 Vibration Equations (Continued) 


= thrust moments about X axis due to engine j 
= thrust moments about Y axis due to engine j 
= thrust moments about Z axis due to engine j 
= slosh moments fcbout X axis at station j 

= slosh moments about Y axis at station j 

= slosh moments about Z axis at station j 

= mode shape translation in X direction for mode i at location j 

= mode shape translation in Y direction for mode i at location j 

= mode shape translation in Z direction for mode i at location j 

= mode slope about X axis for mode i at location j 

-mode slope about Y axis for mode i at location j 
= mode slope about Z axis for mode i at location j 


4.3 Aerodynamic Forces 
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4.3 Aerodynamic Forces (Continued) 


V ■ = Q (Z. - Z ) - R (Y. - Y ) 
axbj w v j eg' ' j eg' 


V aybj = R < X j " X cg> " P < Z j " Z cg> 


V azbj = P ( Y j = V ' Q (X d " V 


», M3 

^akbj = i S . j ^xij q i 


M3 


V .'i = z d> . . q . 
aybj • _ 1 9 yiJ 4 t 


azbj ’ ! ^zij q i 


M3 


IV 


M3 


M3 


axbj s V aybj . l = x *zij q i " V azbj i E = x *yij q i 


IV ^ ' 

V aybj ~ V azbj 1 ^xij q i 


- V 


M3 


axbj i E = 1 * Z 1j 


IV 


M3 


M3 
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4.3 Aerodynamic Forces (Continued) 


v = v * + v " + v 1,1 + v ^ 
v axbj .axbj v axbj v axbj v axbj 


I II III TW 

V . . = V . . + V . . + V . . + V . . 
aybj aybj aybj aybj aybj 


* II 111 T U 

V = V ,. + V ,. + V .. + V . . 
azbj azbj azbj azbj azbj 


The previous 6 equations must be solved simultaneously for V fl 


and v azb j' The solution is as follows: 


let 

a ll 


= -1 
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M3 

Z <t> . . q. 
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M3 
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13 i _ x y yij 


M3 
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zij 
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M3 
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4.3 Aerodynamic Forces (Continued) 


o. = Arctan 

J 


/_ V azbj \ 

( V axbj J 


M. 


a 


h p V 
2 v aj 


F • = Q.S(C . + C .a.) 

axj v XOJ xaj 


F ayj 


q j S C y3j 6 j + C ypj P + q j S C y 6r 6 r 


F • = q . S (C . + C . 6 . ) + q . $ C „ 6 

azj v zoj zaj 4 j zde e 


Where: 




V wxpj * 
V 

wypj’ 

v 

wzpj 


[а] 

[ б ] 

[<$] 


= magnitude of relative wind at Station j 
= azimuth angle of relative wind at Station j 

= relative wind velocity at mass point j in platform 
coordinates 

= transformation described in "coordinate systems" 

= transformation described in "coordinate systems" 

= transformation described in "coordinate systems" 
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4.3 Aerodynamic Forces (Continued) 


V 

V 

V 


axbj 

I 

aybj 

l 

azbj 
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= components of air Velocity in body coordinates 
at. mass point j 


P, Q, R 

V V b 

V V ^cg 


axbj * 

v n 
aybj * 

v * 
V azbj 


W l« 

v axbj * 

\l If, 

aybj * 
v 

“azbj 


/ IV 
axbj’ 

/ IV 
aybj’ 

/ IV 
'azbj 


= angular velocity about X, Y, and Z axes respectively 
= location of mass point j in body coordinates 

= location of center of gravity in body coordinates 


= component of velocity at mass point j due to 
rotation of mass point about c.g. 


= components of velocity of vibrating mass with 
respect to rigid body 


= components of velocity due to perpendicular forces 
being rotated with respect to rigid body due to bending 


M3 




V 

V 

V 


axbj’ 
aybj ’ 
azbj 


= number of bending modes 

= velocity of vehicle relative to the earth in platform 
coordinate 


= components of velocity of mass point j with respect 
to air 
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4.3 Aerodynamic Forces (Continued) 


o . 
J 


M 


V aj 


= angle of attack of mass point j 
= mach number of mass point j 
= velocity of mass point j with respect to air 
= speed of sound 

= dynamic pressure at mass point j 


p 

S 

<5 


mass density of air 
aero reference area 
rudder deflection 
elevator deflection 


C 

xoj’ 

x aj’ 

C 

yej 

C ypj * 

C y«5r’ 
C zoj ’ 
C zaj 5 


A 




= aero coefficients for mass point j 


'z6e / 


= sideslip angle for mass point j 
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4.4 Engine Forces 
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4.4 Engine Forces (Continued) 
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M3 


' . + R + E 4> ' . . q. + 2 c 
YJ v zij 


yj 


. O) 0 1 . + 0' . = 


1=1 Z1J "1 ej ey yj ey yj 


a) 2 e . + y^- (X . - X .) 

ey yj I z2e ecgj ej ; 


e' . = /e' . dt + 0 ' . 
PJ ■ PJ pjo 


e 1 . = /0 . dt + 0 ' . 

PJ PJ pjo 


0' . . = /0 . dt + 0' . 

yj yj yjo 


0' . = /0 . dt + 0 . 

yj yj yjo 


F , = ET. cos e' cos 0' . 

txj j p yj 

F tyj * E tj si " e ;j 

F t 2 j'= - E tj s1 " 6 PJ cc)s 6 yj 




M3 

F exj " ” m ej 

A + 

X 

Z 

i =1 

F eyj " ~ m ej * 

IV 

(X ej 


X1J 


ej ^cg 


'cj eg' 


M3 


< X ej ' X CG> R - < Z ej - Z CG> P + ^ *yij 1l 


- F . . sin e 1 . 

8 xj yj 
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4.4 Engine Forces (Continued) 


F ezj ' " m ej 


r ' .M3 

£ A z " ^ X ej " X cg^ + ^ y ej " P + ? =1 *zij q i 


-F . sin e' • 
exj ° pj 


F txj F txj 

F tyj = F tyj 
F tzj - F izj 


+ F . 
exj 


+ m ej ^ X ej ~ X ecg^ 9 y 


+ m . (X - X .) 0' 
ej 1 ecg ej' p 


+ F 


+ F 


eyj 

ezj 


Where: 


0' . = yaw engine gimbal angle with respect to the mounting 

yj surface for engine j 

e . = yaw engine gimbal angle with respect to rigid body 

coordinates 'for engine j F 



pitch engine gimbal angle with respect to the mounting 
surface for engine j 


e . 
PJ 


*eP 

v>- 

W 


pitch engine gimbal angle with respect to rigid body 
coordinates for engine j 


location of engine j pivot point 
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4.4 Engine Forces (Continued) 

*ecg j ’ 1 

Y_ location of engine j center of mass 

ecgj ■ 

Z ecgj .y 



thrust forces acting on vehicle 


6 pj = er) 9 7 * ne ac t uator srgle 

e . = yaw engine actuator angle 

y j ^ 

e . = pitch engine actuator rate 

r J 

I = moment of inertia of engine bell about Y axis at engine 

gimbal point 

I = moment of inertia of engine bell about Z axis at engine 

gimbal point 

e . = yaw engine actuator rate 

y j 

e pcj = P ltc ^ 9 1 wba 1 angle command 

6 • = yaw gimbal angle command 

ycj 

~ damping factor for pitch engine dynamics 

ep 

S e y = damping factor for yaw engine dynamics 
w, = frequency for engine dynamics 

w a = frequency for actuator dynamics 

0p = present engine rate for pitch 

9^ = present engine rate for yaw 

e = pitch engine angular acceleration 

r 

0^ = yaw engine angular acceleration 
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4.4 Engine Forces (Continued) 


M . = mass of engine j 


A x’ A v’ 

y = linear acceleration of vehicle in body coordinates 
A_ 


X , 

_cg’ 

Y cg , V = vehicle center of gravity 

^cg y 


The characteristicsfrequencies associated with the engine dynamics are 
much higher than the vehicle characteristic frequencies. It is 
recommended that the engine dynamics be integrated separately with an 
integration cycle of 50 milliseconds. 


4.5 Slosh Forces 


M3 • 

* v ,- + 2c e4 ow K* + 0 )%, x v . * - z <p xjj q i - A x - o (Z sJ - Z cg ) 


XJ ^SJ SJ XJ SJ XJ 


i = 1 


+ R (Y . - Y ) 
v sj eg' 


2 _ M3 

■yj **sj “sj "yj J sj yj " z T yu y 

i = 1 


X . + 2c • uj . X . + w*" ■ X . = - <j> . . q . - A - R (X . - X ) 

\/t i ct wt ct -\/i v yi i y ' si Co 7 


+ P (Z sj - Z cg> 
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4.5 Slosh Forces (Continued) 


M3 


X zj + 2 ? sj “sj x zj + 10 sj x zj 


1 ^zij <*1 " A z “ P (Y sj ' ' f cg ) 


+ Q (X c-i " X r „> 

sj eg' 


F sxj " - m sj X xj 


F xyj = “ m sj x yj 


F szj ~ “ m sj x zj 


X xj = /X xj 


x yj = /A yj 


A Zj = A Zj 


X xj " /x xj 


x yi = /x yj 


X zj - A zj 


Where: 

Vj’ x yj’ x 2j 
* • • 
X xj’ x yj’ x zj 

X xj» x yj * x zj 


= displacements of slosh mass j 
= velocities of slosh mass j 
= acceleration of slosh mass j 


236 



4.5 Slosh Forces (Continued) 


C sj 


"sj 


* damping factor for slosh mode j 
= characteristic frequency of slosh mode j 


V 

V SJ’ 

V 


= position of slosh mass j 


«n . = mass of sloshing fluid at mode j 


4.6 Aerodynamic Moments 


^axj F ayj ^ Z cg " Z arj^ " F azj ^ Y cg " Y arj^ 


H ayj ■ F azj < X cg ' x arj> ' F axj < Z cg ' Z arj> 


^azj F axj ^ Y cg " Y arj^ " F ayj ^ X cg " X arj^ 


M axj = S b < C 13J 6 j + C 1 6a aa + C l6r 6r > + * C 1 P p + <h r R) 


n " 


ayj qj s c (C moj + C maj Oj + c « +-c Q c ) 


m<$e e mg 
2 V. 


' b ; 
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4.6 Aerodynami c Moments (Continued) 


M azj = q j S b C ngj + 2 J V axbj C np P + q j S ^ (C n6a 6 a + C n6r 6r) 


M axj M axj + M ax L' 


M ayj M ayj + 


M . + M l + M 
azj azj azj 


Where: 


C l6a’ C n6a’ C lgj’ C l5r* C lpj’ C lr’ C moj’ C jpaj ? . b j’ c j’ W’ C mqj’ C ngj’ 


C„„> ? aero coefficients for station j 

np n<sr ° 


<S 3 = aileron deflection 

a - 


4.7 Engine Moments 


M txj ~ “ T ybj ^ Z ej " Z cg^ + T zbj ^ Y ej " Y cg^ 


M 


W = - T Z bj < X ej - X cg> + T xbj < Z ej ’ Z cg> 


H.- 


* - T xbj < X ej ■ V + T ybj (X ej - V 
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Engine Moments (Continued) 


txj 


It — 


-T, 


M3. M3 

r ybj ^ , \ *zij q 1 + T zbj * , j Vij q i 


"tyj = -T . . M3 » ..q. + T h . « 3 » ..q. 

zbj E 9 xij xbj E 9 zij 


i = 1 


i = 1 


u ii _ T M3 

tzj “xbj E r ylj 
i = 1 


d> . . q . + T M3 ■ 

™ 1 ybj E ^xij q i 


i = 1 


M exj ^ F ezj _ni ej^ X ej~ X ecg^ e *pj ^ ^ Y ej~ Y cg^ F eyj 4tn ej^ X ej” X ecgj 


M eyj fr exj^ Z ej" Z ecgj^^ F ezj’ rn ej^ X ej“ X ecgj^ " e pj^ X ej" X cg^ 


M ezj ^■ F eyj’ Hn ej^ X ej -X ecgj^ e 'yj^ X ej" X cg^" F exj^ Y ej" Y cg^ 


M. . = M ' . + M. '! + M . 

txj txj txj exj 


M tyj = M tyj + M tyj + M eyj 


M tzj = M tzj + M t2 J t M ezj 


Where: 


1 . = moment of inertia of engine j 



4.8 Slosh Moments 


M sxj = A y m sj A zj ' A z A yo 


M syj * A z m sj ^xj ' A x m sj A zi 


M szj = A x W sj A yj ' A y m sj A xi 


M sxj F syj < Z cg “ Z sj> ‘ F sza < Y cg “ V 


M syj * F szj < X cg " X sj> " F sxj < Z cg " V 


M szj = F sxj ^ Y cg “ V ■ F sy^ x cg - V 


M . = M . + M '! 

SXJ S&J SXJ 


M . = M '. + M 
syj syj syj 


M . = M 1 + M 
szj szj szj 
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4.9 


COORDINATE SYSTEMS 


Inertial Pdlar-Equatorial - A right-handed orthogonal system with its origin 
at the center of the earth - X axis in the equatorial plane and positive 
through a reference meridian at the time of liftoff; the reference meri- 
dian is defined by the time of liftoff and the coordinate system used 
for gravity calculations. The Z axis is positive through the North 
Pole. 

Inertial Plumbline - An orthogonal system with its origin at the 
center of the earth, X axis parallel to the launch site gravity vector 
and positive in the direction opposite to gravitational acceleration. 

The Z axis lies in the launch plane and points downrange and the Y 
axis completes a right-handed triad. 

Local Vertical - An orthogonal system with its origin at the center 
of the earth, the X axis, points from the earth center to the vehicle, 
the Z axis is in the plane containing the earth’s rotation axis and 
the X LV axis. The Z axis is perpendicular to the X axis and points 
towards the North Pole. The Y axis completes a right-handed triad. 

Body - An orthogonal system with its origin at the engine gimbal 
pivot plane - X axis positive towards the nose of the vehicle along 
the main propellant tank centerline, Z axis positive "down", and 
the Y axis completes the right-handed system and is positive in the 
direction of the right wing. 
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Transformation matrix from polar-equatorial to plumbline coordinates 



11 

a 12 

a 13 

21 

a 22 

a 23 

31 

a 32 

a 33 


a^ = cos x^ cos {(j) L + w e t L ) 

* / * \ 
a 12 = cos x^ sin (<^> L + <» e t L ) 

-r* 

sin x^ 

•fc "ft 

a 2 j = sin A l sin x L cos (o) e t L +<j> L ) - cos A L sin (w e t|_ + <f>|_) 

•ff ^ 

a 22 = sin A l sin x L sin (u> e t L +<j> L ) - cos A L cos (u> e t L +<j) L ) 

* 

a 23 = -sin A^ cos x^ 

jc *k 

a 31 = _c0 s \ Sln X L cos ^ w e t L +< * > L^ ” sin \ sin ( w e^L +< ^L^ 

1c ic 

a 32 = " cos \ S1 " n X L sin + sin A L cos (w e t L +<{> L ) 

* 


a 33 


cos A l cos x l 


Where: 

**■ 

x^ = geodetic latitude of launch site 

* 

4>^ = longitude of launch site 
w e = angular rate of earth 

= time of launch (from epoch) 

A^ = launch azimuth 
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Transformation matrix from body to inertial plumbline coordinates 


[■]■ 


'll 

b 12 

b 13 

'21 

b 22 

b 23 

'31 

b 32 

b 33 


= cos e cos ip 


'12 


= sin e sin <p - cos e sin \p cos <p 


bis = sin 0 cos 4 > + cos 0 sin tp sin <p 

b 2 j = sin <p 

b 22 = cos <Jj cos <p 


23 


= -cos ip sin <f> 


31 


= -sin 0 cos ip 


- cos 0 sin 4 > + sin 6 sin ip cos <p 


'33 


= cos e cos <p - sin e sin ip sin 


where the Euler angles 0, <P and <j> are calculated in E0M. 



Transformation matrix 


from local vertical to 


polar-equatorial 


coordinates’. 



d ll = 

cos x 

y cos $ 

A - 

-sin 



°12 




d 13 = 

-sin 

X V 

cos $ 

d 21 = 

5 cos 

X V 

sin <j>‘ 

d 22 

= cos 

<J>' 



d 2J * -S 1n X V S ’ n *' 
"31 = Sin 1V 


d 33 = COS x v 


Where: 

sin Xy - ^ — - 

cos X V X V 

sin 4>’ = V R C0S Xv 

cos +' * V R cos X V 
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